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Von Statischen Potenzialen zur Hochenergiestreuung 
Zusammenfassung 



Wir entwickcln cin Loop-Loop Korrclations Modcll zur cinhcitlichen Beschrcibung von statischen 
Farbdipol Potenzialen, einschliessenden QCD Strings, und hadronischen Hochencrgiereaktio- 
nen mit besonderer Berucksichtigung von Sattigungseffekten, die die S-Matrix Unitaritat bei 
extrcm hohen Encrgicn manifestiercn. Das Modell verbindet storungsthcoretischen Gluonaus- 
tausch mit dem nicht-storungstheoretischen Modell des Stochastischen Vakuums, das den Ein- 
schluss von Farbladungen durch Flussschlauchbildung der Farbfelder beschreibt. Wir berech- 
nen Farbfeldverteilungcn statischer Farbdipolc in vcrschiedenen SU(N C ) Darstcllungcn und 
linden Casimir Skalierungsvcrhaltcn in Ubcreinstimmung mit aktuellcn Gitter-QCD Ergeb- 
nissen. Wir untcrsuchen die im einschliessenden String gespeicherte Encrgie und zeigen mit 
Niedcrcnergicthcoremen die Konsistenz mit dem statischen Quark- Antiquark Potential. Wir 
verallgemeinern Meggiolaros Analytische Fortsetzung von Parton-Parton auf Dipol-Dipol Streu- 
ung und erhaltcn cincn Euklidischcn Zugang zur Hochenergiestreuung, der prinzipiell crlaubt, 
Streumatrixelemente in Gitter-QCD zu bcrcchncn. Mit dem Euklidischcn Loop-Loop Korrcla- 
tions Modell berechnen wir in diesen Zugang Dipol-Dipol Streuung bei hohen Encrgicn. Das 
Ergebnis bildet zusammen mit einer universellen Energieabhangigkeit und reaktionsspczifischen 
Wellenfunktionen die Grundlagc fur cine cinheitliche Beschreibung von pp, 7rp, Kp, 7*p und 77 
Reaktionen in guter Ubereinstimmung mit experimentellen Daten fur Wirkungsquerschnittc, 
Steigungsparameter und Strukturfunktionen. Die erhaltenen Stossparamctcrprofile fur pp und 
7£p Reaktionen und die stossparameterabhangige Gluonverteilung des Protons xG{x,Q 2 1 \ b±\) 
zeigen Sattigung bei extrem hohen Energien in Ubereinstimmung mit Unitaritatsgrenzen. 



From Static Potentials to High-Energy Scattering 

Abstract 



We develop a loop-loop correlation model for a unified description of static color dipole po- 
tentials, confining QCD strings, and hadronic high-energy reactions with special emphasis on 
saturation effects manifesting S'-matrix unitarity at ultra-high energies. The model combines 
perturbative gluon exchange with the non-perturbative stochastic vacuum model which describes 
color confinement via flux-tube formation of color fields. We compute the chromo-field distri- 
butions of static color dipoles in various SU(N C ) representations and find Casimir scaling in 
agreement with recent lattice QCD results. We investigate the energy stored in the confining 
string and use low-energy theorems to show consistency with the static quark-antiquark po- 
tential. We generalize Meggiolaro's analytic continuation from parton-parton to dipole-dipolc 
scattering and obtain a Euclidean approach to high-energy scattering that allows us in principle 
to calculate S'-matrix elements in lattice QCD. In this approach we compute high-energy dipolc- 
dipole scattering with the Euclidean loop-loop correlation model. Together with a universal 
energy dependence and reaction-specific wave functions, the result forms the basis for a unified 
description of pp, 7rp, Kp, 7*p, and 77 reactions in good agreement with experimental data 
for cross sections, slope parameters, and structure functions. The obtained impact parameter 
profiles for pp and 7£p reactions and the impact parameter dependent gluon distribution of 
the proton xG(x, Q 2 , \b±\) show saturation at ultra-high energies in accordance with unitarity 
constraints. 
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Chapter 1 



Introduction 



According to our present understanding, quantum chromodynamics (QCD) is the theory 
of strong interactions and thus describes the diversity of strong interaction phenom- 
ena [1]. QCD is the gauge field theory defined by the non-Abelian color gauge group 
SU(N C ) with N c = 3 colors [2] and the presence of a certain number of quark fields ip q 
(q = u,d,s,...) in the fundamental representation. The Lagrange density of QCD has 
an extremely simple form 

£qcd(x) = - \Q%{x)G a ^{x) + J2 U*) {i-fW* + i9Gffi c ) ~ m q ) %{x) (1.1) 

q 

and involves only very few parameters, i.e. the gauge coupling g and the current quark 
masses m q . The first term describes pure gauge theory of the massless gluon potentials 
Q^(x) (a = 1, ...,Nc — 1) in terms of the gluon field strengths 

G%{x) = dnQKx) - dvQ'ix) - gf bc Gl(x)Gt(x) , (1.2) 

where f abc are the structure constants of the gauge group SU(N C ). The second term is 
a sum over the contributions of the quarks with flavor q = u,d,s, ... and involves the 
SU(N C ) generators t" in the fundamental representation r = N c . 

The simple Lagrangian (1.1) brings along a very rich structure. Due to vacuum 
polarization, the effective coupling depends on the distance scale, or equivalently the 
(inverse) energy scale, at which it is measured: As long as there are no more than 
Nf = 16 quark flavors, the renormalization group tells us that the effective coupling 
becomes small at short distances and thus that QCD is an asymptotically free theory [3]. 
Indeed, high-energy deep inelastic scattering experiments reveal that quarks behave as 
free particles at short distances. Accordingly, perturbation theory is applicable in this 
regime and allows reliable analytic calculations, for example, of the total cross section 
for electron-positron annihilation into hadrons [4,5]. At distance scales of order of the 
proton size (« lfm), the effective coupling becomes large so that perturbation theory 
breaks down. 
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Genuinely non-perturbative methods are needed to describe the physics of hadrons 
and low-energy interactions [6]. Indeed, an analytic derivation of color confinement - 
the phenomenon that quarks and gluons cannot be observed as isolated particles - from 
the QCD Lagrangian is still missing and among the ultimate goals in theoretical physics. 
Further fundamental strong interaction phenomena at low energy are spontaneous chiral 
symmetry breaking and dynamical mass generation, i.e. the hadron spectrum and the 
origin of the hadron mass, which are inherently non-perturbative phenomena as well 
not yet proven analytically from the QCD Lagrangian. 

A challenge at high energy is the description and understanding of hadronic high- 
energy scattering [7]. For small momentum transfers, the effective QCD coupling is again 
too large for a reliable perturbative treatment and non-perturbative methods are needed. 
In particular, it is a key issue to unravel the effects of confinement and topologically 
non-trivial gauge field configurations (such as instantons) on such reactions [8,9]. 

The most interesting phenomenon in hadronic high-energy scattering is the rise 
of the total cross sections with increasing cm. energy y/s: While the rise is slow in 
hadronic reactions of large particles such as protons, pions, kaons, or real photons [10], 
it is steep if only one small particle is involved such as an incoming virtual photon [11- 
15] or an outgoing charmonium [16, 17]. This energy behavior is best displayed in the 
proton structure function F2(x,Q 2 ) that is equivalent to the total j*p cross section 
<j*°p(s, Q 2 ). With increasing photon virtuality Q 2 , the increase of F2(x,Q 2 ) towards 
small Bjorken x = Q 2 /s becomes significantly stronger. Together with the steep rise 
of the gluon distribution in the proton xG(x, Q 2 ) with decreasing x, the rise of the 
structure function ^(x, Q 2 ) towards small x [11-15] is one of the most exciting results 
of the HERA experiments. 

In high-energy collisions of stable hadrons, the rise of the total cross sections is 
limited: The Froissart bound, derived from very general principles such as unitarity 
and analyticity of the S-matrix, allows at most a logarithmic energy dependence of 
the cross sections at asymptotic energies [18]. Analogously, the rise of a^ p (s,Q 2 ) is 
expected to slow down. The microscopic picture behind this slow-down is the concept 
of gluon saturation: Since the gluon density in the proton becomes large at high energies 
yfs (small x), gluon fusion processes are expected to tame the growth of <r*** p (s, Q 2 ). It 
is a key issue to determine the energy at which these processes become significant. 

Lattice QCD is the principal theoretical tool to study non-perturbative aspects of 
the QCD Lagrangian from first principles. Numerical simulations of QCD on Euclidean 
lattices give strong evidence for color confinement and spontaneous chiral symmetry 
breaking and describe dynamical mass generation from the QCD Lagrangian [19-21]. 
However, since lattice QCD is limited to the Euclidean formulation of QCD, it cannot be 
applied in Minkowski space-time to simulate high-energy reactions in which particles are 
inherently moving near the light-cone. Furthermore, it is hard to understand from sim- 
ulations the important QCD mechanisms that lead, for example, to color confinement. 
Here (phenomenological) models that allow analytic calculations are important. 
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In this thesis we develop a loop-loop correlation model (LLCM) for a unified de- 
scription of static quark-antiquark potentials and confining QCD strings in Euclidean 
space-time [22] and hadronic high-energy reactions in Minkowski space-time [23, 24] with 
special emphasis on saturation effects manifesting S'-matrix unitarity at ultra-high ener- 
gies [25,26]. The model is constructed in Euclidean space-time. It combines perturba- 
tive gluon exchange with the stochastic vacuum model (SVM) of Dosch and Simonov [27] 
which leads to confinement via a string of color fields [22,28,29]. For applications to 
high-energy scattering, the LLCM can be analytically continued to Minkowski space- 
time [23] following the procedure introduced for applications of the SVM to high-energy 
reactions [30-32]. In this way the LLCM allows us to investigate manifestations of the 
confining QCD string in unintegrated gluon distributions of hadrons and photons [24]. 
We present an alternative Euclidean approach to high-energy scattering that shows 
how one can access high-energy scattering in lattice simulations of QCD. Indeed, this 
approach allows us to compute S'-matrix elements for dipole-dipole scattering in the 
Euclidean LLCM and confirms the analytic continuation of the model to Minkowski 
space-time [22]. The applications of the LLCM to high-energy scattering are based 
on the functional integral approach developed for parton-parton scattering [33, 34] and 
extended to gauge- invariant dipole-dipole scattering [30-32]; see also Chap. 8 of [35]. 
Together with a universal energy dependence introduced phenomenologically, the func- 
tional integral approach to dipole-dipole scattering is the key to the presented unified 
description of hadron-hadron, photon-hadron, and photon-photon reactions. 

The central element in our approach is the gauge-invariant Wegner- Wilson loop [36, 
37]: The considered physical quantities are obtained from the vacuum expectation value 
(VEV) of one Wegner- Wilson loop, (W r [C]), and the correlation of two Wegner- Wilson 
loops, (Wn[Ci]W r2 [C2]). Here ru\ indicates the SU(N C ) representation of the Wegner- 
Wilson loops which we keep as general as possible. In phenomenological applications, 
the propagation of (anti-)quarks requires the fundamental representation, r = N c , and 
the propagation of gluons the adjoint representation, r = N^ — 1. We express (W r [C]) 
and (W ri [Ci]W r2 [C2]) in terms of the gauge-invariant bilocal gluon field strength corre- 
lator integrated over minimal surfaces by using the non-Abelian Stokes theorem and a 
matrix cumulant expansion in the Gaussian approximation. We decompose the gluon 
field strength correlator into a perturbative and a non-perturbative component. Here 
the SVM is used for the non-perturbative low-frequency background field [27] and per- 
turbative gluon exchange for the additional high-frequency contributions. This com- 
bination allows us to describe long and short distance correlations in agreement with 
lattice calculations of the gluon field strength correlator [38,39]. Moreover, it leads to 
a static quark-antiquark potential with color Coulomb behavior for small source sep- 
arations and confining linear rise for large source separations. We calculate the static 
quark-antiquark potential with the LLCM parameters adjusted in fits to high-energy 
scattering data [23] and find good agreement with lattice data. We thus have one model 
that describes both static hadronic properties and high-energy reactions of hadrons and 
photons in good agreement with experimental and lattice QCD data. 
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We apply the LLCM to compute the chromo- electric fields generated by a static color 
dipole in the fundamental and adjoint representation of SU(N C ). The non-perturbative 
SVM component describes the formation of a color flux tube that confines the two 
color sources in the dipole [22, 28, 29] while the perturbative component leads to color 
Coulomb fields. We find Casimir scaling for both the perturbative and non-perturbative 
contributions to the chromo-electric fields in agreement with lattice data and our results 
for the static dipole potential, which is the potential of a static quark-antiquark pair 
in case of the fundamental representation and the potential of a gluino pair in case of 
the adjoint representation. The mean squared radius of the confining QCD string is 
calculated as a function of the dipole size. Transverse and longitudinal energy density 
profiles are provided to study the interplay between perturbative and non-perturbative 
physics for different dipole sizes. The transition from perturbative to string behavior 
is found at source separations of about 0.5 fm in agreement with the recent results of 
Liischer and Weisz [40]. 

The low-energy theorems, known in lattice QCD as Michael sum rules [41], relate the 
energy and action stored in the chromo-fields of a static color dipole to the corresponding 
ground state energy. The Michael sum rules, however, are incomplete in their original 
form [41]. We present the complete energy and action sum rules [42-44] in continuum 
theory taking into account the contributions to the action sum rule found in [29] and the 
trace anomaly contribution to the energy sum rule found in [42] . Using these corrected 
low-energy theorems, we compare the energy and action stored in the confining string 
with the confining part of the static quark-antiquark potential. This allows us to confirm 
consistency of the model results and to determine the values of the (3 - function and the 
strong coupling a s at the renormalization scale at which the non-perturbative SVM 
component is working. Earlier investigations along these lines have been incomplete 
since only the contribution from the traceless part of the energy-momentum tensor has 
been considered in the energy sum rule. 

To study the effect of the confining QCD string examined in Euclidean space-time 
on high-energy reactions in Minkowski space-time, an analytic continuation from Eu- 
clidean to Minkowski space-time is needed. For investigations of high-energy reactions 
in our model constructed in Euclidean space-time, the gauge-invariant bilocal gluon field 
strength correlator can be analytically continued from Euclidean to Minkowski space- 
time. As mentioned above, this analytic continuation has been introduced by Dosch and 
collaborators for applications of the SVM to high-energy reactions [30-32] and is used 
also in our Minkowskian applications of the LLCM [23-26]. Recently, an alternative an- 
alytic continuation for parton-parton scattering has been established in the perturbative 
context by Meggiolaro [45]. This analytic continuation has already been used to access 
high-energy scattering from the supergravity side of the AdS/CFT correspondence [46], 
which requires a positive definite metric in the definition of the minimal surface [47], 
and to examine the effect of instantons in high-energy scattering [48] . 

In this thesis we generalize Meggiolaro's analytic continuation [45] from parton- 
parton to gauge-invariant dipole- dipole scattering such that S'-matrix elements for high- 
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energy reactions can be computed from configurations of Wegner- Wilson loops in Eu- 
clidean space-time and with Euclidean functional integrals. This evidently shows how 
one can access high-energy reactions directly in lattice QCD. First attempts in this 
direction have already been carried out but only very few signals could be extracted, 
while most of the data was dominated by noise [49] . We apply this approach to compute 
the scattering of dipoles at high-energy in the Euclidean LLCM: We recover exactly the 
result derived with the analytic continuation of the gluon field strength correlator [23] . 
This confirms the analytic continuation used in all earlier applications of the SVM 
to high-energy scattering [30-32, 50-61] including the Minkowskian applications of the 
LLCM [23-26]. Here we use the obtained S-matrix element Sdd as the basis for our uni- 
fied description of hadronic high-energy reactions with special emphasis on saturation 
effects in hadronic cross sections and gluon saturation [23,25,26]. Moreover, we have 
used the obtained S'-matrix element to investigate manifestations of the confining string 
in high-energy reactions of hadrons and photons [24]. In particular, we have found that 
the string can be represented as an integral over stringless dipoles with a given dipole 
number density. This decomposition of the confining string into dipoles gives insights 
into the microscopic structure of the model. It allows us to calculate unintegrated gluon 
distributions of hadrons and photons from dipole-hadron and dipole-photon cross sec- 
tions via \k±\ factorization. Our result shows explicitly that non-perturbative physics 
dominates the unintegrated gluon distributions at small transverse momenta \k±\ [24]. 

Aiming at a unified description of hadron-hadron, photon-hadron, and photon- 
photon reactions, we follow the functional integral approach to high-energy scattering 
in the eikonal approximation [30, 31, 33, 34] (cf. also Chap. 8 of [35]) in which S'-matrix 
elements factorize into the universal S-matrix element for elastic high-energy dipole- 
dipole scattering Sod and reaction-specific light-cone wave functions. The color dipoles 
- described by light-like Wegner- Wilson loops - are given by the quark and antiquark 
in the meson or photon and in a simplified picture by a quark and diquark in the 
baryon. Consequently, hadrons and photons are described as color dipoles with size and 
orientation determined by appropriate light-cone wave functions [31, 34]. 

We introduce a phenomenological energy dependence into the univeral S-matrix el- 
ement for dipole-dipole scattering Sdd in order to describe simultaneously the energy 
behavior in hadron-hadron, photon-hadron, and photon-photon reactions involving real 
and virtual photons as well. Motivated by the two-pomeron picture of Donnachie and 
Landshoff [62], we ascribe to our non-perturbative (soft) and perturbative (hard) com- 
ponent a weak and strong energy dependence, respectively. Including multiple gluonic 
interactions, we obtain S-matrix elements with a universal energy dependence that 
respects unitarity constraints in impact parameter space. 

To study saturation effects that manifest S-matrix unitarity, we consider the scat- 
tering amplitudes in impact parameter space, where S-matrix unitarity imposes rigid 
limits on the impact parameter profiles such as the black disc limit. We confirm explic- 
itly that our model respects this unitarity constraint for dipole-dipole scattering, which 
is the underlying process of each considered reaction in our approach. We calculate the 
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impact parameter profiles for proton-proton and longitudinal photon-proton scattering. 
The profile functions describe the blackness or opacity of the interacting particles and 
give an intuitive geometrical picture for the energy dependence of the cross sections. At 
ultra-high energies, the hadron opacity saturates at the black disc limit which tames the 
growth of the hadronic cross sections in agreement with the Froissart bound [18]. We es- 
timate the impact parameter dependent gluon distribution of the proton xG(x, Q 2 , \ b±\) 
from the profile function for longitudinal photon-proton scattering and find gluon satu- 
ration at small Bjorken x that tames the steep rise of the integrated gluon distribution 
xG(x, Q 2 ) towards small x. These saturation effects manifest S'-matrix unitarity in 
hadronic collisions and should be observable in future cosmic ray and accelerator ex- 
periments at ultra-high energies. The cm. energies and Bjorken x at which saturation 
sets in are determined and LHC and THERA predictions are given. 

With the intuitive geometrical picture gained in impact parameter space, we turn 
to experimental observables to analyse the energy dependence of the cross sections and 
to localize saturation effects. We compare the LLCM results with experimental data 
and provide predictions for future cosmic ray and accelerator experiments. Total cross 
sections a tot , the structure function of the proton F2, slope parameters B, differential 
elastic cross sections da el /dt, elastic cross sections a el , and the ratios a el /a tot and 
a tot /B are considered for proton-proton, pion-proton, kaon-proton, photon-proton, and 
photon-photon reactions involving real and virtual photons as well. 

The outline of this thesis is as follows: In Chap. 2 the loop-loop correlation model 
is developed in its Euclidean version and the general computations of (Wj.[C]) and 
(W ri [C\]W r2 [C2]) are presented. Based on these evaluations, we compute in Chap. 3 
potentials and chromo-field distributions of static color dipoles with emphasis on Casimir 
scaling and the interplay between perturbative color Coulomb behavior and non-per- 
turbative formation of the confining QCD string. Moreover, low-energy theorems are 
discussed and used to show consistency of the model results and to determine the values 
of (5 and a s at the renormalization scale at which the non-perturbative SVM component 
is working. In Chap. 4 the Euclidean approach to high-energy scattering is presented 
and applied to compute high-energy dipole-dipole scattering in our Euclidean model. 
The additional ingredients for the unified description of hadronic high-energy scatter- 
ing, i.e. hadron and photon wave functions and the phenomenological universal energy 
dependence, are introduced in Chap. 5. Going to impact parameter space in Chap. 6, 
we confirm the unitarity condition in our model, study the impact parameter profiles 
for proton-proton and photon-proton scattering, and discuss the impact parameter de- 
pendent gluon distribution of the proton xG(x, Q 2 , \ b±\) and gluon saturation. Finally, 
in Chap. 7 we present the phenomenological performance of the LLCM and provide 
predictions for saturation effects in experimental observables. In the Appendices we 
give explicit parametrizations of the loops and the minimal surfaces and provide the 
detailed computations for the results in the main text. 
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Chapter 2 



The Loop-Loop Correlation 



Model 



In this chapter the vacuum expectation value (VEV) of one Wegner- Wilson loop and 
the correlation of two Wegner- Wilson loops are computed for arbitrary loop geome- 
tries within a Gaussian approximation in the gluon field strengths. The results are 
applied in the following chapters. We describe our model for the QCD vacuum in 
which the stochastic vacuum model (SVM) of Dosch and Simonov [27] is used for the 
non-perturbative low-frequency background field (long-distance correlations) and per- 
turbative gluon exchange for the additional high-frequency contributions (short-distance 
correlations) . In this and the next chapter we work in Euclidean space-time as indicated 
by exclusively subscript Dirac indices and space-time variables written in capital letters. 

2.1 Vacuum Expectation Value of a Wegner-Wilson Loop 

A crucial quantity in gauge theories is the Wegner-Wilson loop operator [36, 37] 



Concentrating on SU(N C ) Wegner-Wilson loops, where N c is the number of colors, 
the subscript r indicates a representation of SU(N C ), Tr r (...) = Tr r (...)/Tryll r is the 
normalized trace in the corresponding color space with unit element H r , g is the strong 
coupling, and Gfi(Z) = G^(Z)t^ is the gluon potential with the SU(N C ) group generators 
in the corresponding representation, i", that demand the path ordering indicated by 
V on the closed path C in space-time. A distinguishing theoretical feature of the 
Wegner-Wilson loop is its invariance under local gauge transformations in color space. 
Therefore, it is the basic object in lattice gauge theories [19, 20, 36, 37, 63] and has been 
considered as the fundamental building block for a gauge theory in terms of gauge 




(2.1) 
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invariant variables [64]. Phenomenologically, the Wegner- Wilson loop represents the 
phase factor associated to the propagation of a very massive or very fast color source 
in the representation r of the gauge group SU(N C ). 

To compute the expectation value of a Wegner- Wilson loop (2.1) in the QCD vacuum 



W r [C] 



G 



Tr r V exp 



-igfdZ^Z)^: 



G 



(2.2) 



we transform the line integral over the loop C into an integral over the surface S with 
dS = C by applying the non-Abelian Stokes' theorem [65] 



W r [C] 



G 



Try Vs exp 



-i 9 - da^{Z)G%{0,Z-C zo )t a r 
J s 



G 



(2.3) 



where Vs indicates surface ordering and O is an arbitrary reference point on the sur- 
face S. In Eq. (2.3) the gluon field strength tensor, Q^Z) = Q^ v {Z)t a , is parallel 
transported to the reference point O along the path Czo 



0^(0, Z; Czo) = HO, Z; C Z oY l G^{Z)$>(0, Z- Czo) 
with the QCD Schwinger string 



$(0,Z;C ZO ) = Vexp 



g [ dZ,Q*(Z)t a r 
J Czo 



(2.4) 



(2.5) 



The QCD vacuum expectation value {■ ■ - )g represents functional integrals in which the 
functional integration over the fermion fields has already been carried out as indicated 
by the subscript G [34]. The model we use for the QCD vacuum works in the quenched 
approximation that does not allow string breaking through dynamical quark-antiquark 
production. 



Due to the linearity of the functional integral, (Tr r . . .) = Tr r (. . .), we can write 



W r [C7]) G = Tiv(7>sexp 



i| / da^(Z)g a (0,Z;C zo )t a r 



G 



(2.6) 



For the evaluation of (2.6), a matrix cumulant expansion is used as explained in [34] 
(cf. also [66]) 



Vs exp 



I da(Z)g(0,Z;C ZO ) 
Js 



G 



= ex P[ f T / d*(*i) ■ ■ ■ MXn) K n (X u ■ ■ ■ , X n )} , (2.7) 

n=l 

where space-time indices are suppressed to simplify notation. The cumulants K n consist 
of expectation values of ordered products of the non-commuting matrices 0(O, Z; Czo)- 
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The leading matrix cumulants are 



K X {X) = {G(0,X;C x ))g, 
K 2 (X 1 ,X 2 ) = (P s [Q(0,X 1 ;C Xl )g(0,X 2 ;C X2 )}) G 

~ ((g(0,X 1 ;C Xl )) G (g(0,X 2 ;C X2 )) G + (1 <- 2) 



(2- 



(2.9) 



Since the vacuum does not prefer a specific color direction, K\ vanishes and K 2 becomes 

K 2 {X U X 2 ) = (Vs[g(0,X i; C Xl )g(0,X 2 ;C X2 )}} G . (2.10) 

Now, we approximate the functional integral associated with the expectation values 
(. . .)g as a Gaussian integral in the gluon field strength. Consequently, the cumulants 
factorize into two-point field correlators such that all higher cumulants K n with n > 2 
vanish 1 and (W r [C])G! can be expressed in terms of K 2 



W r [C}) =Tr r exp 

G 



~Y I da^Xx) [ da pa (X 2 ) 
•J S <J s 



s Js 
"Ps[Guv(0, Xi;C Xl o) t a r g b pu {0, X 2 ; C X2 o) t\ 



G. 



(2.11) 



Due to the color neutrality of the vacuum, the gauge-invariant bilocal gluon field 
strength correlator contains a 5-function in color space, 



%v(9i x i;C Xl o)g b pa (0, X 2 ; Cx 2 o) 



='■ F IXI/pa (X 1 ,X 2 ,0;C Xl o,C X2 o) 

(2.12) 

which makes the surface ordering Vs in (2.11) irrelevant. The tensor F pupa will be 
specified in Sec. 2.3. With (2.12) and the quadratic Casimir operator C 2 (r), 



t t 



Eq. (2.11) reads 

(w r [c] 

where 



G 



Tr r exp 



C 2 (r) 



Xss % 



Br , 






[ ° 2{r K 1 


= exp 


2 /voo 



(2.13) 



(2.14) 



Xss J da^X^ J da^(X 2 )F^{X 1 ,X 2 ,0-,C Xl o,C Xa o) ■ (2-15) 

In this rather general result (2.14) obtained directly from the color neutrality of the 
vacuum and the Gaussian approximation in the gluon field strengths, the more detailed 
aspects of the QCD vacuum and the geometry of the considered Wegner- Wilson loop 



1 We are going to use the cumulant expansion in the Gaussian approximation also for perturbative 
gluon exchange. Here certainly the higher cumulants are non-zero. 



9 



are encoded in the function xss which is computed in Appendix B for the rectangular 
loop shown in Fig. 3.1. 

In explicit computations we use for S the minimal surface, which is the planar surface 
spanned by the loop, C = dS, that leads to Wilson's area law [27]. The minimal surface 
is represented in the upcoming figures by the shaded areas (cf. Figs. 3.1 and A.l). Of 
course, the results should not dependent on the surface choice. In our model this will 
be fulfilled for the perturbative and non-perturbative non-confining component but not 
for the non-perturbative confining component in F pupa (specified in Sec. 2.3) due to 
the Gaussian approximation and the associated truncation of the cumulant expansion. 
Nevertheless, since our results for the VEV of a rectangular Wegner-Wilson loop lead 
to a static quark-antiquark potential that is in good agreement with lattice data (see 
Sec. 3.1), we are led to conclude that the choice of the minimal surface is required by 
the Gaussian approximation in the gluon field strengths. The minimal surface is also 
favored by other complementary approaches such as the strong coupling expansion in 
lattice QCD, where plaquettes cover the minimal surface, or large-AT c investigations, 
where the planar gluon diagrams dominate in the large- N c limit. Within bosonic string 
theory, our minimal surface represents the world-sheet of the rigid string: Our model 
does not describe fluctuations or excitations of the string and thus cannot reproduce 
the Liischer term which has recently been confirmed with unprecedented precision by 
Liischer and Weisz [40]. 



2.2 The Loop-Loop Correlation Function 

The computation of the loop-loop correlation function (W ri [C\]W r2 [C 2 ])g starts again 
with the application of the non-Abelian Stokes' theorem [65] that allows us to transform 
the line integrals over the loops Ci j2 into integrals over surfaces S± >2 with dS\p = Ci j2 



ir n [r 1 ]Tr,. 2 [r 2 ]> G = (Tr ri p 5 exp 



i- / da llv (X 1 )g a (0 1 ,X 1 ;C Xl o 1 )tr 1 



Si 



x Tr r2 V s exp 



-i 9 - I da pa (X 2 )g b pa (0 2 ,X 2 ;C X2 o 2 )t 
1 Js 2 



G 



(2.16) 



where 0\ and 2 are the reference points on the surfaces S\ and S 2 , respectively, that 
enter through the non-Abelian Stokes' theorem. In order to ensure gauge invariance in 
our model, the gluon field strengths associated with the loops must be compared at one 
reference point O. Due to this physical constraint, the surfaces S\ and S 2 are required 
to touch at a common reference point 0\ = 2 = O. 

To treat the product of the two traces in (2.16), we transfer the approach of Berger 
and Nachtmann [57], cf. also [23], to Euclidean space-time. Accordingly, the product 
of the two traces respectively over SU(N C ) matrices in the r\ and r 2 representation, 
fr ri (...) Tr r2 (...), is interpreted as one trace Tr ri0r2 (...) := Tr ri0r2 (...)/Tr ri0r2 (lL. 1(g)r . 2 ) 
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that acts in the tensor product space built from the r\ and r 2 representations 



W rx [C x \W r2 \C 2 \ 



G 



TV 



V s exp[-d da llu (XiKu(0,X 1 ;C Xl o)t a ri ] ® % 
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ri 



Psexp[- 



Si 

-i 9 - [da pa (X 2 )g b pa (0,X 2 ;C X2 o)t 
z Js 2 



G 



With the identities 



exp ( t a ri ) (g> ll r2 
ll n <g> exp ( ) 



exp ( ® ]l r2 ) 
exp ( ll n <g> ^ 2 ) 



(2.17) 



(2.18) 
(2.19) 



the tensor products can be shifted into the exponents. Using the matrix multiplication 
relations in the tensor product space 



{t a ri ® H r2 )« ® ll r2 ) = ® ll r2 

.6 



£ ® t; 2 



(2.20) 



® ll r2 ) (ll n ® 4) 
and the vanishing of the commutator 

U,]^®^] =0, (2.21) 
the two exponentials in (2.17) commute and can be written as one exponential 



t a 



(w[d]w[c 2 ] 



Tr ri(glr2 Vs exp 



-if / d^(X)^(0,X;Cxo) 
2 is 



G 



(2.22) 



with the following gluon field strength tensor acting in the tensor product space 



Gnv{0,X;Cxo) 



g« v {0,X;C X o){t a ri ®X 2 ) for X G 5i 
^(O.^Cxo)^! ® C 2 ) for X €. Si 



(2.23) 



In Eq. (2.22) the surface integrals over Si and S2 are written as one integral over the 
combined surface S = S± + S 2 so that the left-hand side (lhs) of (2.22) becomes very 
similar to the lhs of (2.3). This allows us to proceed analogously to the computation 
of (W r [C])G in the previous section. After exploiting the linearity of the functional 
integral, the matrix cumulant expansion is applied, which holds for G(j,v(0,X; Cxo) as 
well. Then, with the color neutrality of the vacuum and by imposing the Gaussian 
approximation now in the color components of the gluon field strength tensor, only the 
n = 2 term of the matrix cumulant expansion survives, which leads to 



W Tl [C x \W r2 \C 2 



= Tr ri0r2 exp 



G 



(2.24) 



7 fda lu/ (X 1 ) [ da pa (X 2 ){Vs[G^(0,X 1 ;C Xl o)G P a(0,X 2 ;Cx 2 o, Jl 
<Js Js /G 
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Note that the Gaussian approximation on the level of the color components of the 
gluon field strength tensor (component factorization) differs from the one on the level 
of the gluon field strength tensor (matrix factorization) used to compute (W r [C]) in the 
original version of the SVM [27]. Nevertheless, with the additional ordering rule [28] 
explained in detail in Sec. 2.4 of [67], a modified component factorization is obtained 
that leads to the same area law as the matrix factorization. 

Using definition (2.23) and relations (2.20), we now redivide the exponent in (2.24) 
into integrals of the ordinary parallel transported gluon field strengths over the separate 
surfaces S\ and S2 



W ri [Cx]W n [C 2 ] > = Tr ri0r2 exp 



(2.25) 

da M 4X 1 )^^ pff (X 2 )P 5 [(^(0,X 1 ;Cx 1 o)^ CT (0,X 2 ;C'x 2 o)) G (^ 1 ® 4 2 ) 

-^Jdv^X^Jdapv^Vs [(g%(0, X i; C Xl0 )Q b pa {0, X 2 ; C X2 oj) G (t a ri t b ri ® % 2 ) 

-£[ da IH ,(X 1 )f da pa (X2)Vs\(g%(0,X 1 ;C Xl0 )g b p AO,X2-,C X2 oj) (ll n ® t a r / r2 ) 
°Js 2 JS 2 LX IG 

Here the surface ordering Vs is again irrelevant due to the color neutrality of the vac- 
uum (2.12), and (2.25) becomes 



W rx [C x \W r2 \C2. 



G 



Tr 



riQ$r 2 



exp 



X.S1S2 X.S 2 S\ 



(t a rl ® t a r2 ) 



XS1S1 (,a + a 



{t a ri t a ri ® % 2 ) 



Xs 2 s 2 



(ll n ® t a r2 t a r2 ) 



(2.26) 



with 



XS iSj / *V(*0 / ^(^2)^(Xi,^2,0;Cx l0 ,Cx 2 o) ■ (2.27) 

*s Si *J Sj 

The symmetries in the tensor structure of F pvpa - see (2.42), (2.44), and (2.48) - lead 
to XS1S2 = XS 2 Si- With the quadratic Casimir operator (2.13) our final Euclidean result 
for general SU (N c ) representations r\ and r 2 becomes 2 



G 



(2.28) 



- XSrS 2 ® - (^pl X5 lSl + Xs 2 s 2 ) Un 



C 2 (r 2 ) 



2 Note that the Euclidean Xs^Si 7^ in contrast to XStSi = for Minkowskian light-like loops Cj 
considered in the original version of the Berger-Nachtmann approach [23,57]. 
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where H ri ®r 2 = \\ ® %-2- After specifying the representations r\ and r 2 , the tensor 
product i n ®r 2 := *ri ® ^r 2 can b e expressed as a sum of projection operators Pj with 
the property Pji ri(X)r2 = AjPj (no sum over i), 



with 



A,; = 



T 1 r ri (g )r2 (Pj in®r 2 ) 



Tr 



ri®r 2 



(2.29) 



which corresponds to the decomposition of the tensor product space into irreducible 
representations. 

For two Wegner- Wilson loops in the fundamental representation of SU(N C ), r\ = 
f2 = N c , that could describe the trajectories of two quark-antiquark pairs, the decom- 
position (2.29) is trivial 



.a 



with the projection operators 



fa _ 



2N r 



-p«- 



N c + 1 
2N r 



Pa , 



(2.30) 



( P s)(aia 2 )(/3i/3 2 ) ~ ^(^ai/3i^a 2 /3 2 + ^ai&^/Ji) 



(Pa)(a 1 a 2 )( / 3 1 /3 2 ) - „ (^ai/Ji <^a 2 /3 2 ~ ^ai/3 2 <J Q2/ 3i ) 



(2.31) 
(2.32) 



that decompose the direct product space of two fundamental SU (N c ) representations 
into the irreducible representations 



N c ® N c = (N c + l)N c /2 N C (N C - l)/2 . 
With TrN c( g)N c Hn c ®n c = N% and the projector properties 



(2.33) 



P 2 Si a = ? s ,a , Tr Nc0Nc P s = (N c + l)JV c /2 and Tr NcS5Nc P a = (N c - l)N c /2 (2.34) 

we find for the loop-loop correlation function with both loops in the fundamental SU(N C ) 
representation 



W Nc [Ci]W Nc [C 2 j/ =exp 

G 



c 2 (n c ; 



XSiSi XS2/S 2 ^ 



(2.35) 



A^c + l 



2N, 



exp 



iVc- 1, 
2N C 



"X5iS 2 



+ 



2iV„ 



exp 



A^c + 1 
2iV n 



'X.S1S2 



where 



N 2 - 1 



(2.36) 



For one Wegner- Wilson loop in the fundamental and one in the adjoint representa- 
tion of SU(N C ), r\ = N c and r 2 = N^ — 1, which is needed in Sec. 3.2 to investigate 
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the chromofield distributions around color sources in the adjoint representation, the 
decomposition (2.29) reads 

*N C ® t a Nhl = - ^ Pi + \ P 2 - \ P 3 (2.37) 

with the projection operators 3 Pi, P2, and P3 that decompose the direct product space 
of one fundamental and one adjoint representation of SU(N C ) into the irreducible rep- 
resentations 

N C ®N2-1 = N c ^N c (N c -l)(N c + 2) © In c (N c + 1)(N c -2) (2.38) 
which reduces for N c = 3 to the well-known SU (3) decomposition 



15 © 6 . 



(2.39) 



With Tr 



N C ®N: 



-1^1 



N c <g>N2-l 



N c (Nc — 1) and projector properties analogous to (2.34), 
we obtain the loop-loop correlation function for one loop in the fundamental and one 
loop in the adjoint representation of SU(N C ) 



W Nc [C l ]W N 2_ 1 [C 2 \) G = CX,> 



- C 2 (Nc) C 2 (N 2 c -l) 

X.SlSl ~T~ X.S2S2 



(2.40) 





\N C 1 


(iv|-i e5Ip 


2 XS1S2 



N c + 2 
+ ^ttt . — exp 



where 



2(iV c + l) 

C 2 (N2-1) = iV c . 



1 

2 XS1S2 



+ 



N c -2 
2(N C -1) 



exp 



2 X5152 



(2.41) 



Note that our expressions for the loop-loop correlation function (2.28) and, more 
specifically, (2.35) and (2.40), are rather general results - as our result for the VEV of one 
Wegner- Wilson loop (2.14) - obtained directly from the color neutrality of the vacuum 
and the Gaussian approximation in the gluon field strengths. The loop geometries, 
which characterize the problem under investigation, are again encoded in the functions 
XSiSj-, where also more detailed aspects of the QCD vacuum enter in terms of F^ upa , 
i.e. the gauge-invariant bilocal gluon field strength correlator (2.12). 

For the explicit computations of XS1S2 presented in Appendix B, one has to specify 
surfaces 5i,2 with the restriction dSip = C\p, according to the non-Abelian Stokes' 
theorem. As illustrated in Figs. 3.4 and 4.2, we choose for S± t 2 the minimal surfaces 
that are built from the planar areas spanned by the corresponding loops C\ )2 and the 
infinitesimally thin tube which connects the two surfaces S\ and S 2 - This is in line 
with our surface choice in applications of the LLCM to high-energy reactions [23-25] 

3 The explicit form of the projection operators Pi, P2, and P3 can be found in [68] but note that we 
use the Gell-Mann (conventional) normalization of the gluons. The eigenvalues, A;, of the projection 
operators in (2.37) can be evaluated conveniently with the computer program "Colour" [69]. 
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illustrated in Fig. 4.1. The thin tube allows us to compare the gluon field strengths in 
surface S\ with the gluon field strengths in surface S^- 

Due to the Gaussian approximation and the associated truncation of the cumu- 
lant expansion, the non-perturbative confining contribution (see Sec. 2.3) to the loop- 
loop correlation function depends on the surface choice. Consequently, our results for 
the chromo-field distributions of color dipoles obtained with the minimal surfaces (see 
Sec. 3.2) differ from the ones obtained with the pyramid mantle choice for the sur- 
faces [28] even if the same parameters are used. With low-energy theorems we show 
in Sec. 3.3 that the minimal surfaces are actually required to ensure the consistency of 
our results for the VEV of one loop, (W r [C]), and the loop- loop correlation function, 

In applications of the model to high-energy scattering [23-26] the surfaces are inter- 
preted as the world-sheets of the confining strings in line with the picture obtained for 
the static dipole potential from the VEV of one loop. The minimal surfaces are the most 
natural choice to examine the scattering of two rigid strings without any fluctuations 
or excitations. Our model does unfortunately not choose the surface dynamically and, 
thus, cannot describe string flips between two non-perturbative color dipoles. Recently, 
new developments towards a dynamical surface choice and a theory for the dynamics of 
the confining strings have been reported [70]. 



2.3 Perturbative and Non-Perturbative QCD Components 

We decompose the gauge-invariant bilocal gluon field strength correlator (2.12) into a 
perturbative (P) and non-perturbative (NP) component 

Ffxvpc = FpLvpo + Ffivpa j (2.42) 

where gives the low-frequency background field contribution modeled by the 

non-perturbative stochastic vacuum model (SVM) [27] and F^ v the additional high- 
frequency contribution described by perturbative gluon exchange. This combination 
allows us to describe long and short distance correlations in agreement with lattice cal- 
culations of the gluon field strength correlator [38,39]. Moreover, this two component 
ansatz leads to the static quark-antiquark potential with color Coulomb behavior for 
small and confining linear rise for large source separations in good agreement with lat- 
tice data as shown in Sec. 3.1. Note that besides our two component ansatz an ongoing 
effort to reconcile the non-perturbative SVM with perturbative gluon exchange has led 
to complementary methods [70-72]. 

We compute the perturbative correlator E% pa from the gluon propagator in Feynman- 
't Hooft gauge 
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where we introduce an effective gluon mass of rxiQ = m p = 0.77 GeV to limit the range of 
the perturbative interaction in the infrared (IR) region. This value is, of course, impor- 
tant for the interplay between the perturbative and non-perturbative component which 
comes out reasonable as illustrated in Sec. 3.1 for the static quark-antiquark poten- 
tial. Moreover, it gives the "perturbative glueball" (GB) generated by our perturbative 
component a reasonable finite mass of M^ B = 2mc = 1.54 GeV. 

In leading order in the strong coupling g, the resulting bilocal gluon field strength 
correlator is gauge-invariant already without the parallel transport to a common refer- 
ence point so that depends only on the difference Z = X± — X 2 



pupa 



(Z) 



g 2 1 



7T" 



d 
dZ v 



\xp 



Z p 5 pa ) + 



d 



dZ ( Zp ^ vc ~ Z °^ v p) D p( Z 



(2.44) 



7T^ 



d 4 K 



iKZ 



(2vr) 4 

with the perturbative correlation function 



K v K a B w - K u K p 5 pa + KpKpSn - K p K a 5 up D' P {K 2 ) 



Dp(Z 2 ) = 
D' P (K 2 ) := 



rn 



G 



2tt 2 Z 2 
d 



K 2 (m,G \ Z\ 



dK 2 



J d 4 Ze lKZ D P (Z 2 ) = - 



K 2 + 



in 



G 



(2.45) 
(2.46) 



The perturbative gluon field strength correlator has also been considered at next-to- 
leading order, where the dependence of the correlator on both the renormalization scale 
and the renormalization scheme becomes explicit and an additional tensor structure 
arises together with a path dependence of the correlator [73]. However, cancellations of 
contributions from this additional tensor structure have been shown [72]. We refer to 
Sec. 3.3 of [67] for a more detailed discussion of this issue. 

We describe the perturbative correlations in our phenomenological applications only 
with the leading tensor structure (2.44) and take into account radiative corrections by 
replacing the constant coupling g 2 with the running coupling 

g 2 (Z 2 ) = 4na s (Z 2 ) = ^ T (2.47) 

(33-2iV»ln [(Z- 2 + M 2 )/A 2 QCD 

in the final step of the computation of the x-function, where the Euclidean distance \Z\ 
over which the correlation occurs provides the renormalization scale. In Eq. (2.47) Nf 
denotes the number of dynamical quark flavors, which is set to Nf = in agreement 
with the quenched approximation, A-qcd = 0.25 GeV, and M allows us to freeze g 2 for 
\Z\ — > 00. Relying on low-energy theorems, we freeze the running coupling at the value 
g 2 = 10.2 (= a s = 0.81), i.e. M = 0.488 GeV, at which our non-perturbative results for 
the confining potential and the total flux tube energy of a static quark-antiquark pair 
coincide (see Sec. 3.3). 
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The tensor structure (2.44) together with the perturbative correlation function (2.45) 
or (2.46) leads to the color Yukawa potential (which reduces for tuq = to the color 
Coulomb potential) as shown in Sec. 3.1. The perturbative contribution thus dominates 
the full potential at small quark-antiquark separations. 

If the path connecting the points X\ and X2 is a straight line, the non-perturbative 
correlator Fj^ depends also only on the difference Z = X\ — X%. Then, the most 
general form of the correlator that respects translational, 0(4), and parity invariance 
reads [27] 

FZAZ) = F^(Z) + F»f™(Z) (2.48) 

= 3(^2 _ -g G 2 (S W 8ua ~ S^Svp) D(Z 2 ) 
1 d d 

- (1 - k) [KvKvS^ - K u K p S^ + KyKpbva - K^K a 5 up ] D[(K 2 )} , 

where 

D[(K 2 ) := J <$ZD^Z 2 )e lKZ . (2.49) 

In all previous applications of the SVM, this form depending only on Z = X\ — X 2 
has been used. New lattice results on the path dependence of the correlator show a 
dominance of the shortest path [74] . This result is effectively incorporated in the model 
since the straight paths dominate in the average over all paths. 

The non-perturbative correlator (2.48) involves the gluon condensate [75] G2 := 

2 

(£^10^(0)0 pu (0)) , the parameter k that determines the non-Abelian character of the 
correlator, and the correlation length a that enters through the non-perturbative corre- 
lation functions D and D\. 

We adopt for our calculations a simple exponential correlation function 

D(Z 2 ) = D X {Z 2 ) = exp(-|Z|/a) , (2.50) 

which is motivated by lattice QCD measurements of the gluon field strength correla- 
tor [38,39]. This correlation function stays positive for all Euclidean distances \Z\ and 
thus is compatible with a spectral representation of the correlation function [76]. This 
means a conceptual improvement since the correlation function that has been used in 
several earlier applications of the SVM [28,29,31,50-58,61] becomes negative at large 
distances. 

With the exponential correlation function (2.50) the lattice data of the gluon field 
strength correlator down to distances of 0.4 fm give the following values for the pa- 
rameters of the non-perturbative correlator [39]: G2 = 0.173 GeV 4 , k = 0.746, and 
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a = 0.219 fm. We have optimized these parameters in our fit to high-energy scattering 
data [23] presented in Chap. 7 (see also Sec. 5.3): 

a = 0.302 fm, k = 0.74, G 2 = 0.074 GeV 4 . (2.51) 

We use these optimized parameters (2.51) throughout this work. They lead to a static 
quark-antiquark potential that is in good agreement with lattice data and, in particular, 
give a QCD string tension of 0-3 = 0.22 GeV 2 = 1.12GeV/fm as shown in Sec. 3.1. 
This value is consistent with hadron spectroscopy [77], Regge theory [78], and lattice 
QCD investigations [63]. Moreover, the non-perturbative component with a = 0.302 fm 
generates a "non-perturbative glueball" with a mass of Mq% = 2/ a = 1.31 GeV which is 
smaller than Mq B = 1.54 GeV and thus governs the long-range correlations as expected. 
We thus have one model that describes both static hadronic properties and high-energy 
reactions of hadrons and photons in good agreement with experimental and lattice QCD 
data. 

Finally, let us emphasize that the non-perturbative correlator (2.48) is a sum of the 
two different tensor structures, F^Jl™ c and F^^£, with characteristic behavior: The 
tensor structure F^™ c is characteristic for Abelian gauge theories, exhibits the same 
tensor structure as the perturbative correlator (2.44) and does not lead to confine- 
ment [27], i.e. it gives an exponentially vanishing static color dipole potential at large 
dipole sizes as shown explictly in Sec. 3.1. In contrast, the tensor structure F^§ can 
only occur in non-Abelian gauge theories and Abelian gauge-theories with monopoles. 
It leads in the case of k 7^ to confinement [27], i.e. to the confining linear increase of 
the static potential at large dipole sizes as demonstrated in Sec. 3.1. Therefore, we call 
the tensor structure multiplied by (1 — k) non-confining (nc) and the one multiplied by 
k confining (c). 
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Chapter 3 



Static Color Dipoles and 
Confining QCD Strings 

In this chapter we apply the loop-loop correlation model to compute the QCD potential 
and the chromo-field distributions of static color dipoles in the fundamental and adjoint 
representation of SU(N C ). Special emphasis is on Casimir scaling behavior and the 
interplay between perturbative Coulomb behavior and non-perturbative formation of 
the confining QCD string. Moreover, low-energy theorems are discussed that relate the 
energy and action stored in the chromo-fields to the static quark-antiquark potential. 
These energy and action sum rules allow us to show consistency of the model results 
and to determine the values of (3 and a s at the renormalization scale at which the 
non-perturbative SVM component is working. 

3.1 The Static Color Dipole Potential 

The static color dipole - two static color sources separated by a distance R in a net 
color singlet state - is described by a Wegner- Wilson loop W r [C] with a rectangular 
path C of spatial extension R and temporal extension T — > oo where r indicates the 
SU(N C ) representation of the considered sources. Figure 3.1 illustrates a static color 
dipole in the fundamental representation r = N c . The potential of the static color dipole 
is obtained from the VEV of the corresponding Wegner- Wilson loop [37, 79] 

V r (R) = -\im Im(WqC]) pot , (3.1) 

1 -^oo 1 

where "pot" indicates the subtraction of the self-energy of the color sources. The static 
quark-antiquark potential Vn c is obtained from a loop in the fundamental representa- 
tion (r = N c ) and the potential of a static gluino pair V^-l from a loop in the adjoint 
representation (r = N^ — 1). 
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With our result for (W r [C]), (2.14), obtained with the Gaussian approximation in 
the gluon field strength, the static potential reads 

V r (R) = ^Q lim ^xsspot, (3.2) 

with the self-energy subtracted, i.e. xsspot '■= Xss — Xss self (see Appendix B). Accord- 
ing to the structure of the gluon field strength correlator, (2.12) and (2.42), there are 
perturbative (P) and non-perturbative (NP) contributions to the static potential 

Vr(R) = ^ T lim I { X p SSpot + (xK + xR) } , (3-3) 

where the explicit form of the x-functions is given in (B.9), (B.28), and (B.37). 

The perturbative contribution to the static potential describes the color Yukawa 
potential (which reduces to the color Coulomb potential [80] for tug = 0) 

V r p (R) = -C 2 (r) ^ exp[-m G i?] . (3.4) 

Here we have used the result for Xss pot gi yen i n (B.37) and the perturbative correlation 
function 

Df\&) := / ^^Wl = -2E!_^M (3.5) 

which is obtained from the massive gluon propagator (2.43). As shown below, the 
perturbative contribution dominates the static potential for small dipoles sizes R. 




Figure 3.1: A static color dipole of size R in the fundamental representation. The rectangular 
path C of spatial extension R and temporal extension T indicates the world-line of the dipole 
described the Wegner- Wilson loop Wn c [C]. The shaded area bounded by the loop C — dS 
represents the minimal surface S used to compute the static dipole potential. 
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The non-perturbative contributions to the static potential, the non-confining com- 
ponent (nc) and the confining component (c), read 



3(iV c 2 - 1) 1 

7T 2 G 2 K f R , 



Vr nc (R) = C 2 (r) ' g^ilx)" &i W {#) (3-6) 



V r NPc (R) = C 2 (r) jf dp(R-p)D^) , (3.7) 

where we have used the results for Xsspot an< ^ XsSpot = Xss ° gi ven respectively 
in (B.28) and (B.9) obtained with the minimal surface, i.e. the planar surface bounded 
by the loop as indicated by the shaded area in Fig. 3.1. With the exponential correlation 
function (2.50), the correlation functions in (3.6) and (3.7) read 

<)(Z 2 ) := J 0^e^D^\K*)5(K 4 ) = - a\Z\* K 2 [\Z\/a] , (3.8) 

D^(Z 2 ) := J^L e tKZ D(K 2 )5(K 4 ) = 2\Z\K 1 [\Z\/a\ . (3.9) 

For large dipole sizes, R > 0.5 fm, the non-confining contribution (3.6) vanishes ex- 
ponentially while the confining contribution (3.7) - as anticipated - leads to confine- 
ment [27], i.e. the confining linear increase, 

V r NP C {R) = a r R + const. . (3.10) 

_R>0.5fm 

Thus, the QCD string tension is given by the confining SVM component [27]: For a 
color dipole in the SU(N C ) representation r, it reads 

o-r = C 2 (r) — ^j— J dZ D{Z ) = C 2 (r) — , (3.11) 

where the exponential correlation function (2.50) is used in the final step. Since the 
string tension can be computed from first principles within lattice QCD [63], rela- 
tion (3.11) puts an important constraint on the three parameters of the non-perturbative 
QCD vacuum a, G 2 , and k. With the values for a, G 2 , and k given in (2.51), that are 
used throughout this work, one obtains for the string tension of the SU(3) quark- 
antiquark potential (r = 3) a reasonable value of 

o- 3 = 0.22 GeV 2 = 1.12 GeV/fm . (3.12) 



The static SU(N C = 3) quark-antiquark potential Vn c (R) = V^lR) is shown as a 
function of the quark-antiquark separation R in Fig. 3.2, where the solid, dotted, and 
dashed lines indicate the full static potential and its perturbative and non-perturbative 
contributions, respectively. For small quark-antiquark separations R < 0.5 fm, the per- 
turbative contribution dominates giving rise to the well-known color Coulomb behavior. 
For medium and large quark-antiquark separations R > 0.5 fm, the non-perturbative 
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Figure 3.2: The static SU(N C = 3) quark-antiquark potential V^ c (R) — V 3 (R) as a function 
of the quark-antiquark separation R. The solid, dotted, and dashed lines indicate the full static 
potential and its perturbative and non-perturbative contributions, respectively. For small quark- 
antiquark separations, R < 0.5 fm, the perturbative contribution dominates and gives rise to the 
well-known color Coulomb behavior at small distances. For medium and large quark-antiquark 
separations, R > 0.5 fm, the non-perturbative contribution dominates and leads to the confining 
linear rise of the static potential. As our model is working in the quenched approximation, string 
breaking cannot be described, which is expected to stop the linear increase for R > 1 f m [63, 81]. 

contribution dominates and leads to the confining linear rise of the static potential. The 
transition from perturbative to string behavior takes place at source separations of about 
0.5 fm in agreement with the recent results of Liischer and Weisz [40]. This supports 
our value for the gluon mass mg = m p = 0.77 GeV which is only important around 
R ~ 0.4 fm, i.e. for the interplay between perturbative and non-perturbative physics. 
For R < 0.3 fm and R > 0.5 fm, the effect of the gluon mass, introduced as an IR regu- 
lator in our perturbative component, is negligible. String breaking is expected to stop 
the linear increase for R > 1 fm where lattice investigations show deviations from the 
linear rise in full QCD [63, 81]. As our model is working in the quenched approximation, 
string breaking through dynamical quark-antiquark production is excluded. 

As can be seen from (3.2), the static potential shows Casimir scaling which emerges 
in our approach as a trivial consequence of the Gaussian approximation used to truncate 
the cumulant expansion (2.7). Indeed, the Casimir scaling hypothesis [82] has been 
verified to high accuracy for SU(3) on the lattice [83,84] (see also Fig. 3.3). These 
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Figure 3.3: The static SU(N C — 3) potential of color dipoles in the fundamental representation 
V^i?) (solid line) and adjoint representation Vg(ii) (dashed line) as a function of the dipole size 
R in comparison to ST/ (3) lattice data for (3 = 6.0, 6.2, and 6.4 [63, 84]. The model results are in 
good agreement with the lattice data. This particularly holds for the obtained Casimir scaling 
behavior. 



lattice results have been interpreted as a strong hint towards Gaussian dominance in 
the QCD vacuum and thus as evidence for a strong suppression of higher cumulant 
contributions [85, 86]. In contrast to our model, the instanton model can neither describe 
Casimir scaling [86] nor the linear rise of the confining potential [87]. 

Figure 3.3 shows the static SU(N C = 3) potential for fundamental sources Vn c (-R) = 
Vs(R) (solid line) and adjoint sources V^^R) = Vs(R) (dashed line) as a function of 
the dipole size R in comparison to 51/(3) lattice data [63,84]. The model results are 
in good agreement with the lattice data. In particular, the obtained Casimir scaling 
behavior is strongly supported by SU(3) lattice data [83,84]. This, however, points 
also to a shortcoming of our model: From Eq. (3.2) and Fig. 3.3 it is clear that string 
breaking is neither described for fundamental nor for adjoint dipoles in our model which 
indicates that not only dynamical fermions (quenched approximation) are missing but 
also some gluon dynamics. 
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3.2 Chromo-Field Distributions of Color Dipoles 



As already explained in Sec. 3.1, the static color dipole - two static color sources sepa- 
rated by a distance R in a net color singlet state - is described by a Wegner- Wilson loop 
W r [C] with a rectangular path C of spatial extension R and temporal extension T — ► oo 
(cf. Fig. 3.1) where r indicates the SU(N C ) representation of the considered sources. A 
second small quadratic loop or plaquette in the fundamental representation placed at 
the space-time point X with side length Rp — > and oriented along the af3-ax.es 



Pf(X)=Ti Nc e W 



-ig f dZ^{Z)t% c 



i-R%^g a aP {x)g^{x) + o{R%) 



(3.13) 

is needed - as a "Hall probe" - to calculate the chromo-field distributions at the space- 
time point X caused by the static sources [88, 89] 



AGf. a/3 (X) -.= (^^g^ j3 {x)g^ j3 {x) 



W r [C] 



v 1 N c 
= — nm — -. — » 

Rp^O R% IT 2 



^POC,PO) vac (3-14) 

(3.15) 



(Wr[C\Pg(X)) „ 



with no summation over a and (3 in (3.13), (3.14), and (3.15). In definition (3.14) 
(■ ■ -)w r [C] indicates the VEV in the presence of the static color dipole while (. . .) vac 
indicates the VEV in the absence of any color sources. Depending on the plaquette 
orientation indicated by a and (3, one obtains from (3.15) the squared components of 
the chromo-electric and chromo-magnetic field at the space-time point X 



AG 2 ra p(X) - 



( o 







B y 
B 2 



B 2 X 
E v Bl 



El\ 

El 
J 



(X) , 



(3.16) 



i.e. space-time plaquettes (a(3 = iA) measure chromo-electric fields and space-space 
plaquettes (af3 = ij) chromo-magnetic fields. As shown in Fig. 3.4, we place the static 
color sources on the Ai-axis at (X\ = ±i?/2, 0, 0, X4) and use the following notation 
plausible from symmetry arguments 



E 2 = Ez 



Ei 



El 



Bl = Bi 



Bl 



(3.17) 



Figure 3.4 illustrates also the plaquette P^f(X) at X = (Xi,X 2 , 0, 0) needed to compute 
E?AX). Due to symmetry arguments, the complete information on the chromo-field 
distributions is obtained from plaquettes in "transverse" space X = (X±, X2,0, 0) with 
four different orientations, af3 = 14, 24, 13, 23, cf. (3.17). 
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Figure 3.4: The plaquette-loop geometry needed to compute the squared chromo-electric field 
E?AX) generated by a static color dipole in the fundamental SU(N C ) representation (r = N c ). 
The rectangular path C indicates the world-line of the static dipole described the Wegner- Wilson 
loop WnJC]. The square with side length Rp illustrates the plaquette P^ 4 (X). The shaded 
areas represent the minimal surfaces used in our computation of the chromo-field distributions. 
The thin tube allows us to compare the gluon field strengths in surface Sp with the gluon field 
strengths in surface Sw- 



The energy and action density distributions around a static color dipole in the 
SU(N C ) representation r are given by the squared chromo-field distributions 



with signs according to Euclidean space-time conventions. Low-energy theorems that 
relate the energy and action stored in the chromo-fields of the static color dipole to the 
corresponding ground state energy are discussed in the next section. 

For the chromo-field distributions of a static color dipole in the fundamental rep- 
resentation of SU(N C ), i.e. a static quark-antiquark pair, we obtain with our results 
for the VEV of one loop (2.14) and the correlation of two loops in the fundamental 




(3.18) 



(3.19) 
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representation (2.35) 



AG 2 (X) 



1 N c 
hm — T — exp 

Rp^o R% it 2 



G 2 (N C 



' "KSpSp 



N c + 1 

2N r 



■ exp 



2N r 



~"X.Sp Sw 



+ 



N c - 1 

2N r 



exp 



2N r 



(3.20) 



- 1 



where xSiS, is denned in (2.27). The subscripts P and W indicate surface integrations 
to be performed over the surfaces spanned by the plaquette and the Wegner- Wilson- 
loop, respectively. Choosing the surfaces - as illustrated by the shaded areas in Fig. 3.4 
- to be the minimal surfaces connected by an infinitesimal thin tube (which gives no 
contribution to the integrals) it is clear that xs P s P °c Rp and Xs P s w °c Rp. Being 
interested in the chromo-fields at the space-time point X, the extension of the quadratic 
plaquette is taken to be innnitesimaHy small, Rp — > 0, so that one can expand the 
exponential functions and keep only the term of lowest order in Rp 



AG 2 (X) = -C 2 (N C 



lim 1 1 v 2 

Rp™0 ~R\ 47T 2 " F 



(3.21) 



This result - obtained with the matrix cumulant expansion in a very straightforward 
way - agrees exactly with the result derived in [28] with the expansion method. Indeed, 
the expansion method agrees for small ^-functions with the matrix cumulant expansion 
(Berger-Nachtmann approach) used in this work but breaks down for large x-functions, 
where the matrix cumulant expansion is still applicable. 

The chromo-field distributions of a static color dipole in the adjoint representation 
of SU(N C ), i.e. a static gluino pair, are computed analogously. Using our result (2.40) 
for the correlation of one loop in the fundamental representation (plaquette) with one 
loop in the adjoint representation (static sources), one obtains 



AG N2-1 a/3 



C 2 (N C 



~X.Sp Sp 



1 



exp 



N r 



' XSpSw 



N c +2 



1 



' XSpSw 



+ 



N c -2 
2(JV C -1) 



exp 



1 



XSpSw 



(3.22) 



which reduces - as explained for sources in the fundamental representation - to 

AGL_ la JX) = -C 2 (Nl-l) Jim ^-^x! P 5 w . (3.23) 



Thus, the squared chromo-electric fields of an adjoint dipole differ from those of a fun- 
damental dipole only in the eigenvalue of the corresponding quadratic Casimir operator 
C 2 {r). In fact, Casimir scaling of the chromo-field distributions holds for dipoles in any 
representation r of SU (N c ) in our model. As can be seen with the low-energy theorems 
discussed below, this is in line with the Casimir scaling of the static dipole potential 
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found in the previous section. Besides lattice investigations that show Casimir scal- 
ing of the static dipole potential to high accuracy in SU(3) [83, 84], Casimir scaling of 
the chromo-field distributions has been considered on the lattice as well but only for 
SU{2) [90]. Here only slight deviations from the Casimir scaling hypothesis have been 
found that were interpreted as hints towards adjoint quark screening. 

In our model the shape of the field distributions around the color dipole is iden- 
tical for all SU(N C ) representations r and given by Xs P s w - ^ n * s a S am illustrates the 
shortcoming of our model discussed in the previous section. Working in the quenched 
approximation, one expects a difference between fundamental and adjoint dipoles: string 
breaking cannot occur in fundamental dipoles as dynamical quark-antiquark production 
is excluded but should be present for adjoint dipoles because of gluonic vacuum polariza- 
tion. Comparing (3.21) with (3.23) it is clear that this difference is not described in our 
model. In fact, as shown in Sec. 3.1, string breaking is neither described for fundamen- 
tal nor for adjoint dipoles. Interestingly, even on the lattice there has been no striking 
evidence for adjoint quark screening in quenched QCD [91]. It is even conjectured that 
the Wegner- Wilson loop operator is not suited to studies of string breaking [92]. 

In the LLCM there are perturbative (P) and non-perturbative (NP) contributions to 
the chromo-electric fields according to the structure of the gluon field strength correlator, 
(2.12) and (2.42), 



{XSpS w ( X )) af3 + {XSpSw ( X )) a /3 + UspS w ( X )) a/3 

where we have demanded the non-interference of perturbative and non-perturbative 
correlations in line with the Minkowskian applications of our model [23-26]. In the 
following we give only the final results of the x - functions for the minimal surfaces 
shown in Fig. 3.4. Details on their derivation can be found in Appendix B. 

The perturbative contribution (P) described by massive gluon exchange leads, of 
course, to the well-known color Yukawa field that reduces to the color Coulomb field 
for m g = 0. It contributes only to the chromo-electric fields, E 2 = E 2 (a/3 = 14) and 
E\ = El = E 2 Z [aj3 = 24), and reads explicitly for X =(X 1 ,X 2 , 0, 0) 

c>2 poo 

(x P SpSw (X)) u = --£J dr{(X 1 -R/2)g 2 {Zl A )D P {Zl A ) 

- (X 1 + R/2) g\z\ c ) Dp(Z 2 lc )} (3.25) 

c>2 poo 

(x P SpSw (X)) M = --f-J dTX 2 {g\zl A )D P (Zl A )-g\Zl c )D P {Zt c )} (3.26) 
with the perturbative correlation function (2.45), the running coupling (2.47), and 

Z\ A =(x 1 -?pj +X 2 2 + t 2 and Z 2 C = (x x + f ) + A + r 2 . (3.27) 
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The non-confining non-perturbative contribution (NPnc) has the same structure 
as the perturbative contribution - as expected from the identical tensor structure - 
but differs, of course, in the prefactors and the correlation function, D\ / D p . Its 
contributions to the chromo-electric fields E 2 = E 2 (aj3 = 14) and E\ = E 2 = E 2 
{a/3 = 24) read for X = (X 1 ,X 2 , 0, 0) 

(x»w) M = - %w-i) K) jy{ {xi " m di{z1a) 

- (X 1 + R/2) D 1 (Z 2 C )} (3.28) 

u& c (*)) 24 = - ^M=^ Sj tX ^ Di{z1a) - Di{zlc) } (3 - 29) 

with the exponential correlation function (2.50) and Z 2 A and Z\ c as given in (3.27). 

The confining non-perturbative contribution (NP c) has a different structure that 
leads to confinement and flux-tube formation. It gives only contributions to the chromo- 
electric field E 2 = El {a(3 = 14) that read for X = (X 1 ,X 2 , 0, 0) 

(X»(*)) M = ^^Tyf^^^ 2 )' ( 3 - 3 °) 

with the correlation function given in (3.9) as derived from the exponential correlation 
function (2.50), and 

Z\ = [X x + (1/2 - p)R} 2 + X\ . (3.31) 

In our model there are no contributions to the chromo-magnetic fields, i.e. the static 
color charges do not affect the magnetic background field 

= B 2 X = and B\ = B 2 y = B 2 Z = , (3.32) 

which can be seen from the corresponding plaquette-loop geometries as pointed out in 
Appendix B. Thus, the energy and action densities are identical in our approach and 
completely determined by the squared chromo-electric fields 

e r (X) = s r (X) = -\e 2 {X) . (3.33) 

This picture is in agreement with other effective theories of confinement such as the 
't Hooft-Mandelstam picture [93] or dual QCD [94] and, indeed, a relation between the 
dual Abelian Higgs model and the SVM has been established [95]. In contrast, lattice 
investigations work at scales at which the chromo-electric and chromo-magnetic fields 
are of similar magnitude [44,96]. Using low-energy theorems, we will see in the next 
section, that the vanishing of the chromo-magnetic fields determines the value of the 
/3-function at the renormalization scale at which the non-perturbative component of our 
model is working. 
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Figure 3.5: Energy density distributions g 2 ez(Xi, X 2 — X3) generated by a color dipole in the 
fundamental SU(3) representation (r = 3) for quark-antiquark separations of R = 0.1, 0.5, 1 and 
1.5 fm. Flux-tube formation leads to the confining QCD string with increasing dipole size R. 

In Fig. 3.5 the energy density distributions g 2 e^(Xi, X2 = X^) generated by a color 
dipole in the fundamental SU (3) representation (r = 3) are shown for quark-antiquark 
separations of R = 0.1, 0.5, 1 and 1.5 fm. With increasing dipole size R, one sees 
explicitly the formation of the flux tube which represents the confining QCD string. 

The longitudinal and transverse energy density profiles generated by a color dipole 
in the fundamental representation (r = 3) of SU(N C = 3) are shown for quark-antiquark 
separations (dipole sizes) of R = 0.1, 0.5, 1 and 1.5 fm in Figs. 3.6 and 3.7. The per- 
turbative and non-perturbative contributions are given in the dotted and dashed lines, 
respectively, and the sum of both in the solid lines. The open and filled circles indicate 
the quark and antiquark positions. As can be seen from (3.15) and (3.16), we cannot 
compute the energy density separately but only the product g 2 e r {X). Nevertheless, 
a comparison of the total energy stored in chromo-electric fields to the ground state 
energy of the color dipole via low-energy theorems yields g 2 = 10.2 (= a s = 0.81) for 
the non-perturbative SVM component as shown in the next section. 

In Figs. 3.6 and 3.7 the formation of the confining string (flux tube) with increasing 
source separations R can again be seen explicitly: For small dipoles, R = 0.1 fm, per- 
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Figure 3.6: Longitudinal energy density profiles g 2 £z(Xi, = X$ = 0) generated by a color 
dipole in the fundamental SU(3) representation (r = 3) for quark-antiquark separations of R = 
0.1, 0.5, 1 and 1.5 fm. The dotted and dashed lines give the perturbative and non-perturbative 
contributions, respectively, and the solid lines the sum of both. The open and filled circles 
indicate the quark and antiquark positions. For small dipoles, R = 0.1 fm, perturbative physics 
dominates and non-perturbative correlations are negligible. For large dipoles, R > 1 fm, the 
formation the confining string (flux tube) can be seen which dominates the chromo-electric 
fields between the color sources. 



30 



20 



> 
II 



15 
10 



15 



10 - 



I 1 1 1 1 I 

Total 

Pcrturbativc 
Non-Pcrturbativc 




\ 1 ' ' ' r 

Total 

Pcrturbativc 
Non-Pcrturbativc 



Total 

Pcrturbativc 
Non-Perturbative 



H = 1.0fin 




H=1.5fm 




Figure 3.7: Transverse energy density profiles g 2 es(X 2 , X\ — X 3 — 0) generated by a color 
dipole in the fundamental SU(3) representation (r = 3) for quark-antiquark separations of R = 
0.1, 0.5, 1 and 1.5 fm. The dotted and dashed lines give the perturbative and non-perturbative 
contributions, respectively, and the solid lines the sum of both. The filled circles indicate the 
positions of the color sources. For small dipoles, R — 0.1 fm, perturbative physics dominates 
and non-perturbative correlations are negligible. For large dipoles, R > 1 fm, the formation the 
confining string (flux tube) can be seen which dominates the chromo-electric fields between the 
color sources. 
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Figure 3.8: Root mean squared radius R ms of the flux tube and energy density in the center 
of a fundamental SU(3) dipole g 2 e-i(X — 0) as a function of the dipole size R. Perturbative 
and non-perturbative contributions are given respectively in the dotted and dashed lines and 
the sum of both in the solid lines. For large R, both the width and height of the flux tube in the 
central region are governed completely by non-perturbative physics and saturate respectively at 
-Rm7°° ~ 0.55 fm and e^°°{X = 0) ps 1 GeV/fm 3 . The latter value is extracted with the result 
g 2 = 10.2 deduced from low-energy theorems in the next section. 



turbative physics dominates and non-perturbative correlations are negligible. For large 
dipoles, R > 1 fm, the non-perturbative correlations lead to formation of a narrow flux 
tube which dominates the chromo-electric fields between the color sources. 

Figure 3.8 shows the evolution of the transverse width (upper plot) and height (lower 
plot) of the flux tube in the central region of the Wegner- Wilson loop as a function of 
the dipole size R where perturbative and non-perturbative contributions are given in 
the dotted and dashed lines, respectively, and the sum of both in the solid lines. The 
width of the flux tube is best described by the root mean squared (ms) radius 



JdX ± Xlg*e r (X 1 =0,X ± ) 
Hms ~ V / dX ± X ± g*e r (X x = 0, X ± ) ' ^ 

which is universal for dipoles in all SU(N C ) representations r as the Casimir factors 
divide out. The height of the flux tube is given by the energy density in the center of 
the considered dipole, g 2 e r (X = 0). For large source separations, R > lfm, both the 
width and height of the flux tube in the central region of the Wegner- Wilson loop are 
governed completely by non-perturbative physics and saturate for a fundamental SU (3) 
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dipole (r = N c = 3) at reasonable values of 

^m7°° ~ 0.55 fm and e§^°°(X = 0) « 1 GeV/fm 3 with <? 2 = 10.2 . (3.35) 

Note that the qualitative features of the non-perturbative SVM component do not 
depend on the specific choice for the parameters, surfaces, and correlation functions 
and have already been discussed with the pyramid mantle choice of the surface and 
different correlation functions in the first investigation of flux-tube formation in the 
SVM [28]. The quantitative results, however, are sensitive to the parameter values, the 
surface choice, and the correlation functions and are presented above with the LLCM 
parameters, the minimal surfaces, and the exponential correlation function [22]. 



3.3 Low-Energy Theorems 

Many low-energy theorems have been derived in continuum theory by Novikov, Shifman, 
Vainshtein, and Zakharov [97] and in lattice gauge theory by Michael [41]. Here we 
consider the energy and action sum rules - known in lattice QCD as Michael sum rules 
- that relate the energy and action stored in the chromo-fields of a static color dipole 
to the corresponding ground state energy [37, 79] 

E r (R) = - lim -\n(W r [C]} . (3.36) 

In their original form [41], however, the Michael sum rules are incomplete [29,42]. In 
particular, significant contributions to the energy sum rule from the trace anomaly of the 
energy-momentum tensor have been found [42] that modify the naively expected relation 
in line with the importance of the trace anomaly found for hadron masses [98] . Taking all 
these contributions into account, the energy and action sum rule read respectively [42- 
44] 

E r {R) = j d 3 Xe r (X) - j d 3 Xs r (X) (3.37) 

E r (R)+R d -^ = -^M J d'Xs r (X) (3.38) 

where (3(g) = {idgjd[i with the renormalization scale [i. Inserting (3.38) into (3.37), we 
find the following relation between the total energy stored in the chromo-fields £y ot (i?) 
and the ground state energy E r (R) 



Er 0t (R) := J d 3 Xe r (X) = ^ (^3E r {R) - R 



dEJR)\ 



The difference from the naive expectation that the full ground state energy of the static 
color sources is stored in the chromo-fields is due to the trace anomaly contribution [42] 
described by the second term on the right-hand side (rhs) of (3.37). 
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With the low energy theorems (3.38) and (3.39) the ratio of the integrated squared 
chromo-magnetic to the integrated squared chromo-electric field distributions can be 
derived 

. = !#XB\X) = (2 + 6 0(g) /g) E r {R) + (1 - (5(g) /g) R ggg 
fd?X&(X) (2 -6 /%)/</) E r (R) + (1 + 0(g)/g) R 8 -^ ' 

which becomes for E r (R) = a r R + E se \f after subtraction of the self-energy contribu- 
tions, i.e. the linear potential V r (R) = a r R with string tension a r in the considered 
representation r, 

Q(R) = l + P }fl / / 9 ■ (3-41) 

VV ; v r (R)=* r R 2- (3(g)/ 'g K ' 

In our model there are no contributions to the chromo-magnetic fields (3.32) so that - 
as already discussed in the previous section - the energy and action densities are identical 
and completely determined by the squared chromo-electric fields (3.33). Since the non- 
perturbative SVM component of our model describes the confining linear potential for 
large source separations R, this allows us to determine from (3.41) immediately the 
value of the (3 - function at the scale \x np at which the non-perturbative component is 
working 

M 



9 



= -2 . (3.42) 



Concentrating on the confining non-perturbative component (NPc) we now use 
(3.39) to determine the value of a s = g 2 /(An) at which the non-perturbative SVM com- 
ponent is working. The rhs of (3.39) is obtained directly from the confining contribution 
to the static potential E? Pc (R) = V r NPc (R) given in (3.7). The lhs of (3.39), however, 
involves a division by the a priori unknown value of g 2 after integrating g 2 e r (X) for the 
chromo-electric field of the confining non-perturbative component (3.30). As discussed 
in the previous section, we cannot compute the energy density separately but only the 
product g 2 e r (X). Adjusting the value of g 2 such that (3.39) is exactly fulfilled for source 
separations of R = 1.5 fm, we find that the non-perturbative component is working at 
the scale fi^p at which 

g 2 (HNp) = 10.2 = a s (fi NP ) = 0.81 . (3.43) 

As already mentioned in Sec. 2.3, we use this value as a practical asymptotic limit for the 
simple one-loop coupling (2.47) used in our perturbative component. Note that earlier 
SVM investigations along these lines have found a smaller value of a s (fiNp) = 0.57 with 
the pyramid mantle choice for the surface [28, 29] but were incomplete since only the 
contribution from the traceless part of the energy-momentum tensor has been considered 
in the energy sum rule. 

In Fig. 3.9 we show the total energy stored in the chromo-field distributions around 
a static color dipole in the fundamental (r = 3) and adjoint (r = 8) representation of 
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Figure 3.9: The total energy stored in the chromo-field distributions around a static color 
dipolc of size R in the fundamental (r — 3) and adjoint (r = 8) representation of SU(3) 
from the confining non- perturbative SVM component, °{R), for a s = 0.81 (solid lines) 

compared with the relation to the corresponding ground state energy (dashed lines) given by 
the low-energy theorem (3.39). Good consistency is found even down to very small values of R. 

SU(3) from the confining non- perturbative SVM component, El°^ NPc (R), for a s = 0.81 
(solid lines) as a function of the dipole size R. Comparing this total energy, which 
appears on the lhs of (3.39), with the corresponding rhs of (3.39) (dashed lines), we 
find good consistency even down to very small values of R. This is a nontrivial and 
important result as it confirms the consistency of our loop-loop correlation result - 
needed to compute the chromo-electric field - with the result obtained for the VEV of 
one loop - needed to compute the static potential V r NPc (R). Moreover, it shows that 
the minimal surfaces ensure the consistency of our non-perturbative component. The 
good consistency found for the pyramid mantle choice of the surface relies on the naively 
expected energy sum rule [28, 29] in which the contribution from the traceless part of 
the energy- momentum tensor is not taken into account. 
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Chapter 4 



Euclidean Approach to 
High-Energy Scattering 



In this chapter we present a Euclidean approach to high-energy reactions of color dipoles 
in the eikonal approximation. We give a short review of the functional integral approach 
to high-energy scattering, which is the basis for the presented Euclidean approach and 
for our investigations of hadronic high-energy reactions in the following chapters. We 
generalize the analytic continuation introduced by Meggiolaro [45] from parton-parton 
scattering to dipole-dipole scattering. This shows how one can access high-energy re- 
actions directly in lattice QCD. We apply this approach to compute the scattering of 
dipoles in the fundamental and adjoint representation of SU (N c ) at high-energy in the 
Euclidean LLCM. The result shows the consistency with the analytic continuation of 
the gluon field strength correlator used in all earlier applications of the SVM and LLCM 
to high-energy scattering. Finally, we comment on the QCD van der Waals potential 
which appears in the limiting case of two static color dipoles. 



4.1 Functional Integral Approach to High-Energy Scattering 

In Minkowski space-time high-energy reactions of color dipoles in the eikonal approxi- 
mation are considered - as basis for hadron-hadron, photon-hadron, and photon-photon 
reactions - in the functional integral approach to high-energy collisions developed origi- 
nally for parton-parton scattering [33, 34] and then extended to gauge-invariant dipole- 
dipole scattering [30-32]. The corresponding T-matrix element for the elastic scattering 
of two color dipoles at transverse momentum transfer q\ (t = — cfj 2 ) and cm. energy 
squared s reads 

T r M 2 (s,t,z 1 ,r 1± ,z 2 ,?2±) = 2is JdH^^ [l - S^ r2 (s,b ± , Zl ,r 1± ,z 2 ,r 2± )\ (4.1) 
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with the S-matrix element (M refers to Minkowski space-time) 



5^ r2 (s,6 ± ,zi,r 1± ,z 2 ,r 2± ) = lim x + L 1 2V " . (4.2) 

T^oo (l4/ ri [GiJ)M\l / » / r 2 [L'2j;M 

The color dipoles are considered in the SU(N C ) representation rj and have transverse 
size and orientation r*jj_. The longitudinal momentum fraction carried by the quark of 
dipole i is Z{ . (Here and in the following we use several times the term quark generically 
for color sources in arbitrary SU(N C ) representations.) The impact parameter between 
the dipoles is [51] 

b± = r lq + (1 - zi)r lL - r 2q ~ (1 - z 2 )f 2 ± = r lcm - r 2cm , (4.3) 

where ?i q (fig) is the transverse position of the quark (antiquark), fi± = fi q — fi q , and 
^icm = Zifi q + (1 — Zi)fig is the center of light-cone momenta. Figure 4.1 illustrates the 
(a) space-time and (b) transverse arrangement of the dipoles. The dipole trajectories Cj 
defining the Wegner- Wilson loops in (4.2) are described as straight lines. This is a good 
approximation as long as the kinematical assumption behind the eikonal approximation, 
s » — t, holds that allows us to neglect the change of the dipole velocities Vi = pi/m 
in the scattering process, where pi is the momentum and m the mass of the considered 
dipole. Moreover, the paths Ci are considered light-like 1 in line with the high-energy 
limit, m? <C s — > oo. For the hyperbolic angle or rapidity gap between the dipole 
trajectories 7 = {v\ ■ v 2 ) - which is the central quantity in the analytic continuation 
discussed below and also defined through s = 4m 2 cosh 2 (7/2) - the high-energy limit 
implies 

lim 7 f» ln(s/m 2 ) — > 00 . (4.4) 

The QCD VEVs (. . .)m in the S'-matrix element (4.2) represent Minkowskian functional 
integrals [34] in which - as in the Euclidean case discussed above - the functional 
integration over the fermion fields has already been carried out. 

The S'-matrix element Sp D := S^ Nc for the scattering of light-like dipoles in the 
fundamental SU (N c = 3) representation (77 = r 2 = N c = 3) is the key to our unified 
description of hadron-hadron, photon-hadron, and photon-photon reactions in the fol- 
lowing chapters. With color dipoles given by the quark and antiquark in the meson or 
photon or in a simplified picture by a quark and diquark in the baryon, we describe 
hadrons and photons as quark-antiquark or quark-diquark systems, i.e. fundamental 
SU(3) dipoles, with size and orientation determined by appropriate light-cone wave 
functions [31, 32]. Accordingly, the T-matrix element for the reaction ab — > cal factorizes 
into the universal S-matrix element S|f D and reaction-specific light-cone wave functions 



1 In fact, exactly light-like trajectories (7 — » 00) are considered in most applications of the functional 
integral approach to high-energy collisions [23-26, 30-32, 50-57, 59-61]. A detailed investigation of the 
more general case of finite rapidity 7 can be found in [58] . 
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Figure 4.1: High-energy dipole-dipole scattering in the eikonal approximation represented by 
Wegner- Wilson loops in the fundamental representation of SU(N C ): (a) Space-time and (b) 
transverse arrangement of the Wegner- Wilson loops. The shaded areas represent the strings 
extending from the quark to the antiquark path in each color dipole. The thin tube allows us 
to compare the gluon field strengths in surface Si with the gluon field strengths in surface S2. 
The impact parameter connects the centers of light-cone momenta of the dipoles. 
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V> a ,ft and i^c^d that describe the fj and Zi distribution of the color dipoles [31, 32, 34] 

T ab ^ cd (s,t) = 2is Jd 2 b ± e 1 ^ Jd Zl d 2 ri Jdz 2 d 2 r 2 (4.5) 

x ipt{zi,fi±)ij;* t (z2,f2±) 1- Sn D (s,b ± ,z 1 ,r 1± ,z 2 ,r 2 ±) ip a {z 1 , f 1± ) ip b (z 2 , f 2 ±) . 

Concentrating in this work on reactions with a = c and b = d, the squared wave functions 
|-0iOi,n_L)| 2 := ipc( z iin±) ipa(zi,r 1± ) and \ip 2 (z 2 , r 2 j_)| 2 := ip d (z 2 , r 2 ±) tpb(z 2 , r 2i _) arc 
needed. We use for hadrons the phenomenological Gaussian wave function [61, 106] 
and for photons the perturbatively derived wave function with running quark masses 
m,f(Q 2 ) to account for the non-perturbative region of low photon virtuality Q 2 [52]. In 
Sec. 5.1 we specify and discuss these wave functions explicitly. The scattering of dipoles 
with fixed size \ri\ and fixed longitudinal quark momentum fraction Zi averaged over all 
orientations, 

I D i | 

is considered in Sec. 6.1 to show that 5-matrix unitarity constraints are respected in our 
model. For the analytic continuation of high-energy scattering to Euclidean space-time, 
we now return to the scattering of dipoles with fixed size and orientation r*j and fixed 
longitudinal quark momentum fraction z%. 



4.2 Analytic Continuation of Dipole-Dipole Scattering 

The Euclidean approach to the described elastic scattering of dipoles in the eikonal 
approximation is based on Meggiolaro 's analytic continuation of the high-energy parton- 
parton scattering amplitude [45]. Meggiolaro's analytic continuation has been derived in 
the functional integral approach to high-energy collisions [33, 34] in which parton-parton 
scattering is described in terms of Wegner- Wilson lines: The Minkowskian amplitude, 
(/ M (7,T, t), given by the expectation value of two Wegner- Wilson lines, forming an 
hyperbolic angle 7 in Minkowski space-time, and the Euclidean "amplitude," g E (Q, T, t), 
given by the expectation value of two Wegner- Wilson lines, forming an angle 6 6 [0, n] 
in Euclidean space-time, are connected by the following analytic continuation in the 
angular variables and the temporal extension T, which is needed as an IR regulator in 
the case of Wegner- Wilson lines, 

g E (@,T,t) = g M ( 7 ^iQ : T^-iT,t) , (4.7) 
5 M ( 7 ,T,t) = gE(e^-i>y,T^iT,t) . (4.8) 

Generalizing this relation to gauge-invariant dipole-dipole scattering described in terms 
of Wegner- Wilson loops [31, 32, 34], the IR divergence known from the case of Wegner- 
Wilson lines vanishes and no finite IR regulator T is necessary. Thus, the Minkowskian 
5-matrix element (4.2), given by the expectation values of two Wegner- Wilson loops, 
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forming an hyperbolic angle 7 in Minkowski space-time, can be computed from the 
Euclidean "S-matrix element" 



S rir2 (6, b±, zi , ri_i_, z 2 , f 2 ±) 



r^(W ri [C 1 ]) E (W r2 [C 2 ]) E 



lim 



(4.9) 



given by the expectation values of two Wegner- Wilson loops, forming an angle G £ [0, ir] 
in Euclidean space-time, via an analytic continuation in the angular variable 

S ^r 2 {l ~ ^[s/m 2 ],b±,z 1 ,r 1± ,z 2 ,f 2± ) = Sr ir2 (@ -> -i^f,b ± ,z 1 ,r 1 _ L ,z 2 ,r 2 _ L ) , (4.10) 

where E indicates Euclidean space-time and the QCD VEVs (. . .)e represent Euclidean 
functional integrals that are equivalent to the ones denoted by {■ ■ -)g i n the preceding 
sections, i.e. in which the functional integration over the fermion fields has already been 
carried out. 

The angle is best illustrated in the relation of the Euclidean S'-matrix element (4.9) 
to the van der Waals potential between two static dipoles V rir2 (Q = 0, 6, zi, ri, z 2 , f 2 ), 
discussed in Sec. 4.4, 



Sf r (0, b± , zi , fi± , z 2 , f 2 ± ) 



lim exp 



TV rir2 (Q,b±,z 1 ,f 1 ±,z 2 ,f 2 ±) 



(4.11) 



Figure 4.2 shows the loop-loop geometry necessary to compute S^ r2 (6 / 0, ■ ■ • ) and how 
it is obtained by generalizing the geometry relevant for the computation of the potential 
between two static dipoles (0 = 0): While the potential between two static dipoles is 
computed from two loops along parallel "temporal" unit vectors, t\ = t 2 = (0,0,0, 1), 
the Euclidean S'-matrix element (4.9) involves the tilting of one of the two loops, e.g. 
the tilting of t\ by the angle towards the X3-axis, t\ = (0, 0, — sin 0, cos 0). The 
"temporal" unit vectors U are also discussed in Appendix A together with another 
illustration of the tilting angle 0. 

Since the Euclidean S-matrix element (4.9) involves only configurations of Wegner- 
Wilson loops in Euclidean space-time and Euclidean functional integrals, it can be com- 
puted directly on a Euclidean lattice. First attempts in this direction have been carried 
out but only very few signals could be extracted, while most of the data was dominated 
by noise [49]. Once precise results are available, the analytic continuation (4.10) will 
allow us to access hadronic high-energy reactions directly in lattice QCD, i.e. within a 
non-perturbative description of QCD from first principles. More generally, the presented 
gauge-invariant analytic continuation (4.10) makes any approach limited to a Euclidean 
formulation of the theory applicable for investigations of high-energy reactions. Indeed, 
Meggiorlaro's analytic continuation has already been used to access high-energy scat- 
tering from the supergravity side of the AdS / CFT correspondence [46] , which requires 
a positive definite metric in the definition of the minimal surface [47], and to examine 
the effect of instantons on high-energy scattering [48]. 
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Figure 4.2: The loop-loop geometry necessary to compute 5^ r2 (0 ^ 0, • • • ) illustrated as a 
generalization of the geometry relevant for the computation of the van der Waals potential 
between two static dipolcs (0 = 0). While the potential between two static dipoles is computed 
from two loops along parallel "temporal" unit vectors, t\ = t% = (0,0,0,1), the Euclidean S- 
matrix element (4.9) involves the tilting of one of the two loops, e.g. the tilting of t\ by the 
angle towards the X 3 -axis, t\ = (0, 0, — sin©, cos 0). 



4.3 Dipole-Dipole Scattering in the LLCM 



Let us now perform the analytic continuation explicitly in our Euclidean loop-loop cor- 
relation model. For the scattering of two color dipoles in the fundamental representation 
of SU(N C ), the Euclidean S'-matrix element becomes with the VEVs (2.14) and (2.35) 



S^ D {Q,bx,zi,ri±,Z2,r 2 ±) 



lim 



N c +1 

2N r 



■ exp 



5 'n c n c ( > h,zi,f 1± , z 2 ,f 2± ) 



N c - 1 

2N r 



-XS!S 2 



+ 



N c -1 

2N r 



exp 



2N r 



(4.12) 



where xSiSj ~ defined in (2.27) - decomposes into a perturbative (P) and non-per- 
turbative (NP) component according to our decomposition of the gluon field strength 
correlator (2.42), 



XS1S2 — XS1S1 + XS1S2 — XS1S2 + \Xs-lS2 + XS1S2) ■ 
In the limit T\ = T 2 = T — > 00 and for € [0, 7r], the components read 



x£iS 2 = cot x p 



NP nc 
XS-1S2 



cot x 



NP nc 



x^ = cote x NPc 



(4.13) 



(4.14) 
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with 



X 



X 



NP nc 



X 



NP c 



g'D'W { \f lq _ r - 2 -|) + g* D f) _ ? 2q \) 
- g 2 Df ] (\f lq - r 2q \) - g 2 D'P (\r lq - r 2q 

7T 2 G 2 (1 - K 

3(iV c 2 " 1) 

TT 2 G 2 K 



D[ {2) (\r lq - r 2 ,\) + D'^ (\r lq - r 2q \) 
- D'^ (\r lq - r 2q \) - D'^ (\f lq - r 2q \) 



(4.15) 



(4.16) 



3(iV 2 - 1) 



(fi • r 2 ) 



/ dvi / dv 2 D (2) (\ri q +v 1 f 1 ±-r 2q -v 2 r2±\) (4.17) 
Jo Jo 



as derived explicitly in Appendix B with the minimal surfaces illustrated in Fig. 4.2. 
In Eq. (4.15) the shorthand notation g 2 D'j!s 2 \\Z ± \) = g 2 (\Z ± \) Dp (2) (|Zl|) is used with 
9 2 (\Z±\) again understood as the running coupling (2.47). The transverse Euclidean 
correlation functions 



D^{Z 2 ) := / 



(2vr) 2 



e tKZ D x (K 2 )5(K 3 )5(K A ) 



(4.18) 



are obtained from the (massive) gluon propagator (2.43) and the exponential correlation 
function (2.50) 



Df\z 2 



: ±-K (m G \Z ± \) 

Df\zl) = vra 4 (3+3^ + ^ 
V a a z 

rt 2 \Z 2 ) = 2vra 2 (l + ^)exp 



exp 



\Z±\ 



_]*!) 

a / 



(4.19) 
(4.20) 
(4.21) 



With the full 6-dependence exposed in (4.14), the analytic continuation (4.10) reads 

X5iS 2 = cot6x g -> -il ; cot(-i7)x g -> oo , «x (4.22) 

and leads to the desired Minkowskian S'-matrix element for elastic dipole-dipole scat- 
tering (DD) in the high-energy limit in which the dipoles move on the light-cone 



lim Sff D (s,b±,zi ,r 1± ,z 2 ,r 2 ±) := lim S¥ N (s, b±, Zi, f 1± , z 2 , r 2 ±) 



5^ (cot© -> i,b±,z 1 ,r 1± ,z 2 ,r 2± ) 



lim 



N c + 1 



■ exp 



.N c -1 



+ 



A^c-1 
2N r 



■ exp 



; iVc+ l 



-x 



(4.23) 



where x = X P + nc + X NF c with (4.15), (4.16), and (4.17). 



.NP 
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It is striking that exactly the same result has been obtained in [23] 2 with the alter- 
native analytic continuation introduced for applications of the SVM to high-energy re- 
actions [30-32]. In this complementary approach the gauge-invariant bilocal gluon field 
strength correlator is analytically continued from Euclidean to Minkowskian space-time 
by the substitution <5 MP — > — g^ p and the analytic continuation of the Euclidean corre- 
lation functions to real time D^(Z 2 ) — > D^{z 2 ). In the subsequent steps, one finds 
(W[C])m = 1 due to the light-likeness of the loops and that the longitudinal correlations 
can be integrated out (W ri [C\]W r2 [C2])m = f{s,b±, • • • ). One is left with exactly the 
Euclidean correlations in transverse space that have been obtained above. This con- 
firms the analytic continuation used in the earlier LLCM investigations in Minkowski 
space-time [23-26] and in all earlier SVM applications to high-energy scattering [30- 
32,50-58]. 

In the limit of small ^-functions, \x P \ ^ 1 an d Ix^l ^ 1) (4.23) reduces to 

hm S% D {s, b ± , Zl , r 1± , z 2 , r 2± ) « 1 + ^± = 1 + ■ ( 424 ) 

The perturbative correlations, (x P ) 2 , describe the well-known two-gluon exchange con- 
tribution [99, 100] to dipole-dipole scattering, which is, of course, an important suc- 
cessful cross-check of the presented Euclidean approach to high-energy scattering. The 
non-perturbative correlations, (x NP ) 2 , describe the corresponding non-perturbative two- 
point interactions that contain contributions of the confining QCD strings to dipole- 
dipole scattering. We have analyzed these string contributions systematically as man- 
ifestations of confinement in high-energy scattering reactions and have indeed found 
a new characteristic structure (different from the perturbative dipole factors) in mo- 
mentum space [24]. This analysis has also shown explicitly that the non-perturbative 
contribution governs - as expected - the region of low transverse momenta \k±\. Here, 
we focus on the structure in space-time representation and refer the reader for comple- 
mentary insights to our momentum-space analysis [24]. 

As evident from the v\ and v 2 integrations in (4.17) and Fig. 4.1b, there are con- 
tributions from the transverse projections of the minimal surfaces (Si, 2X1. connecting 
the quark and antiquark in each of the two dipoles. These are the contributions that 
we interpret as manifestations of the strings confining the quark and antiquark in each 
dipole. We thus understand the confining component Xc P as a string-string interaction. 
Interestingly, we have found in dipole-hadron and dipole-photon interactions that the 
strings confining the quark-antiquark pair in the dipole can represented as an integral 
over stringless dipoles with a given dipole number density. As already mentioned, this 
decomposition of the confining string into dipoles even allows us to compute uninte- 
grated gluon distributions of hadrons and photons and thus gives new insights into the 
microscopic structure of the non-perturbative SVM [24]. 

2 To see this identity, recall that (VF[C]) = 1 for light-like loops and consider in [23] the result (2.30) 
for the loop-loop correlation function (2.3) together with the x-function (2.40) and its components given 
in (2.49), (2.54), and (2.57) with the transverse Minkowskian correlation functions (2.50), (2.55), and 
(2.58). 
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Both non-perturbative components, Xc P an d Xnc \ show color transparency for small 
dipoles, i.e. a dipole-dipole cross section with ctdd (^1,^2) oc \ri\ 2 \f 2 \ 2 for |ri2| — ► 0, as 
known for the perturbative case [101]. This can be seen by squaring (4.16) and (4.17) 
to obtain the leading terms in the T-matrix element for small dipoles; see (4.24). 

Due to the truncation of the cumulant expansion in the Gaussian approximation, 
a considerable dependence of Xc P on the specific surface choice is observed. In fact, a 
different and more complicated result for Xc P was obtained with the pyramid mantle 
choice for the surfaces Si 2 in earlier applications of the SVM to high-energy scatter- 
ing [31,32,50-57,61]. However, we use minimal surfaces in line with our model ap- 
plications in Euclidean space-time discussed in the previous chapter. Moreover, the 
simplicity of the minimal surfaces allows us to give an analytic expression for the lead- 
ing term of the non-perturbative dipole-dipole cross section [24]. Phenomenologically, 
in comparison with pyramid mantles, the description of the slope parameter B(s), the 
differential elastic cross section da el /dt(s,t), and the elastic cross section a el {s) can 
be improved with minimal surfaces as shown in Chap. 7. In contrast to the confin- 
ing component Xc^ '> ^ ne n on-confining components, Xnc ano - X P ■> depend only on the 
transverse position between the quark and antiquark of the two dipoles and are therefore 
independent of the surface choice. 

With the insights from the small- x limit, one sees clearly that the full S-matrix 
element (4.23) describes multiple gluonic interactions. Indeed, the higher order terms 
in the expansion of the exponential functions in (4.23) are crucial to respect S-matrix 
unitarity constraints in impact parameter space [23, 57] as shown explicitly in Chap. 6. 

Concerning the energy dependence, the S-matrix element (4.23) leads to energy- 
independent cross sections in contradiction to the experimental observation. Although 
disappointing from the phenomenological point of view, this is not surprising since 
our approach does not describe explicit gluon radiation needed for a non-trivial en- 
ergy dependence. However, based on the S-matrix element (4.23), a phenomenological 
energy dependence can be constructed - see Sec. 5.2 - that allows a unified descrip- 
tion of high-energy hadron-hadron, photon-hadron, and photon-photon reactions and 
an investigation of saturation effects in hadronic cross sections manifesting S-matrix 
unitarity [23, 25, 26]. This, of course, can only be an intermediate step. For a more fun- 
damental understanding of hadronic high-energy reactions in our model, gluon radiation 
and quantum evolution have to be implemented explicitly. 

Note that x = Xc P + Xnc + X P is a real-valued function. Since, in addition, the 
wave functions \ipi{zi, n_i_)| 2 used in this work - see Sec. 5.1 - are invariant under the 
replacement (fi± — > — fi±,Zi — > 1 — Zi), the T-matrix element (4.5) with Sjj D given 
in (4.23) becomes purely imaginary and reads for N c = 3 

T(s,t) = 2is Jd% ± e^ Jdz 1 d\ 1 Jdz 2 d\ 2 \Mzi,ri±)\ 2 \Mz2,r 2 ±)\ 2 (4.25) 
1 - ^ cos (^x(b±,zi,f 1± , z 2 , f 2 ±)\ - ^ cos f ^X{b±,z 1 ,f 1± ,z 2 ,f 2± ) 
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The real part averages out in the integration over r*j and z% because of 
X{b±, 1 - zi, -r 1± ,z 2 ,f 2 ±) = ~x(b±, zi,r 1± , z 2 ,r 2 ±) , 



(4.26) 



which can be seen directly from (4. 15), (4. 16), and (4.17) as (r*ij_ — > — ?i±,zi — ► 1 — zi) 
implies ri g — > rig. In physical terms, (r*jj_ — > — fi±,Zi — > 1 — Zj) corresponds to charge 
conjugation, i.e. the replacement of each parton with its antiparton and the associated 
reversal of the loop direction. Consequently, the T-matrix (4.25) describes only charge 
conjugation C = +1 exchange. Since in our quenched approximation purely gluonic 
interactions are modeled, (4.25) describes only pomeron 3 but not reggeon exchange. 

Although the scattering of two color dipoles in the fundamental representation of 
SU (N c ) is the most relevant case, we can derive immediately also the Minkowskian S- 
matrix element for the scattering of a fundamental (D) and an adjoint dipole ("glueball" 
GB) in the Euclidean LLCM. Using (2.40) and proceeding otherwise as above, we find 
in the high-energy limit 



lim S DGB (s,b, zi,fi±,z 2 ,f 2 ±) := lim 5 N n2 , (9, b, z 1: fi±, z 2 , f 2 ±) (4.27) 

s— >oo S^OO "c'V x 



lim 



h- 1 


r.Jv c 1 




I Y 

. 2 



N c +2 



l 2 X 



+ 



N c -2 
2{N c -l 



■ exp 



i - Y 
2 



where X = X P + X NPnc + X NPc with (4.15), (4.16), and (4.17). 



4.4 Comments on the QCD van der Waals Potential 

Finally, we would like to comment on the QCD van der Waals interaction between two 
color dipoles, which is - as mentioned together with (4.11) - related to the Euclidean 
S'-matrix element in the limiting case of = 0: The QCD van der Waals potential 
between two static dipoles reads in terms of Wegner- Wilson loops [102, 103] 

"-<e = 0, i, „ = 1/2, ft, „ = 1/2, ft) = - & I b • (4.28) 

In this limit (G = 0) intermediate octet states and their limited lifetime become im- 
portant as is well known from perturbative computations of the QCD van der Waals 
potential between two static color dipoles [102-104]: Working with static dipoles, i.e. 
infinitely heavy color sources, there is an energy degeneracy between the intermediate 
octet states and the initial (final) singlet states that leads for perturbative two-gluon 

3 Oddoron C = — 1 exchange is excluded in our model. It would survive in the following cases: (a) 
Wave functions are used that are not invariant under the transformation (f» — » — — » 1 — Zi). (b) 
The proton is described as a system of three quarks with finite separations modeled by three loops 
with one common light-like line, (c) The Gaussian approximation that enforces the truncation of the 
cumulant expansion is relaxed and additional higher cumulants are taken into account. 
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exchange to a linear divergence in T as T — > oo. This IR divergence can be lifted by 
introducing manually an energy gap between the singlet ground state and the excited 
octet state and thus a limit on the lifetime of the intermediate octet state [102-104]. 

In the perturbative limit of g 2 — > and T large but finite, i.e. \ P ^ 1) the perturba- 
tive component of our model describes the two-gluon exchange contribution to the van 
der Waals potential which is plagued by the discussed IR divergence resulting from the 
static limit. In the more general case of g 2 finite and T — > oo, which does not exclude 
non-perturbative physics, one cannot use the small- x limit and multiple gluonic inter- 
actions become important. Here our perturbative component describes multiple gluon 
exchanges that reduce to an effective one-gluon exchange contribution to the van der 
Waals potential whose interaction range (oc l/rag) contradicts the common expecta- 
tions. Indeed, it is also in contradiction to our results for the glueball mass Mqb which 
determines the interaction range (oc 1/Mqb) between two color dipoles for large dipole 
separations. As already mentioned in Sec. 2.3, we find for the perturbative component, 
Mq B = 2rriG, i.e. half of the interaction range of one-gluon exchange, by computing the 
exponential decay of the correlation of two small quadratic loops Prf for large Euclidean 
times r — ► oo 

M OB : = -l. m iln <*f<°>*y» . (4.29) 
t-oot « S (0)){P,f (T)> 

Note that we find for the non-perturbative component, Mq B = 2/ a, which is smaller 
than Mq B = 2rriG with the LLCM parameters and thus governs the long range correla- 
tions in the LLCM. 

Thus, for a meaningful investigation of the QCD van der Waals forces within our 
model, one has to go beyond the static limit in order to describe the limited lifetime of 
the intermediate octet states appropriately. This we postpone for future work since the 
focus in this work is on high-energy scattering where the gluons are always exchanged 
within a short time interval due to the light-likeness of the scattered particles and the 
finite correlation lengths. Nevertheless, going beyond the static limit in the dipole-dipolc 
potential means going beyond the eikonal approximation in high-energy scattering and 
it is, of course, of utmost importance to see how such generalizations alter our results. 
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Chapter 5 



Hadronic Wave Functions and 
Universal Energy Dependence 



In this chapter hadron and photon wave functions are provided and a universal energy 
dependence is constructed. Together with the T-matrix element (4.25), these ingredients 
are crucial for our unified description of hadron-hadron, photon-hadron, and photon- 
photon reactions in the following chapters. The model parameters adjusted in fits to 
experimental data are summarized at the end of this chapter. 

5.1 Hadron and Photon Wave Functions 

The light-cone wave functions ipi(zi,fi) provide the distribution of transverse size and 
orientation r*j and longitudinal quark momentum fraction z; L to the light-like Wegner- 
Wilson loops W[Cj] that represent the scattering color dipoles. In this way, they specify 
the projectiles as mesons, baryons described as quark-diquark systems, or photons that 
fluctuate into a quark-antiquark pair before the interaction. 

The Hadron Wave Function 

In this work mesons and baryons are assumed to have a quark-antiquark and quark- 
diquark valence structure, respectively. As quark-diquark systems are equivalent to 
quark-antiquark systems [105], this allows us to model not only mesons but also baryons 
as color dipoles represented by Wegner- Wilson loops. To characterize mesons and 
baryons, we use the phenomenological Gaussian Wirbel-Stech-Bauer ansatz [106] 




(5.1) 
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Table 5.1: Hadron Parameters 



Hadron 




S h [fm] 


p,p 


0.3 


0.86 




2 


0.607 


K± 


0.57 


0.55 



where the hadron wave function normalization to unity 

Jdzirfn \rl>i(zi,ri)\ 2 = 1 , (5.2) 

requires the normalization constant 

N h = ^ dz t z t (l - zi) e -fe-|) 2 /(2A^) (5 3) 

Jo 

The different hadrons considered - protons, pions, and kaons - are specified by Azh 
and Sh providing the width for the distributions of the longitudinal momentum fraction 
carried by the quark Zi and transverse spatial extension \fi\, respectively. In this work 
the extension parameter Sh is a fit parameter that should resemble approximately the 
electromagnetic radius of the corresponding hadron [61], while Azh = w/(y/2mh) [106] 
is fixed by the hadron mass rrih and the value w = 0.35 — 0.5 GeV extracted from 
experimental data. We find for (anti-)protons Az p = 0.3 and S p = 0.86 fm, for pions 
Az w = 2 and S n = 0.607 fm, and for kaons Azk = 0.57 and Sk = 0.55 fm which are the 
values used in the main text. For convenience they are summarized in Table 5.1. 

Concerning the quark-diquark structure of the baryons, the more conventional three- 
quark structure of a baryon would complicate the model significantly but would lead 
to similar predictions once the model parameters are readjusted [31]. In fact, there 
are also physical arguments that favor the quark-diquark structure of the baryon such 
as the SI = 1/2 enhancement in semi-leptonic decays of baryons [105] and the strong 
attraction in the scalar diquark channel in the instanton vacuum [107]. 

The Photon Wave Function 

The photon wave function ip~/(zi, fi, Q 2 ) describes the fluctuation of a photon with vir- 
tually Q 2 into a quark-antiquark pair with longitudinal quark momentum fraction z% 
and spatial transverse size and orientation f\. The computation of the corresponding 
transition amplitude {qq(zi, Ti)\j*(Q 2 )) can be performed conveniently in light-cone per- 
turbation theory [108] and leads to the following squared wave functions for transverse 
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(T) and longitudinally (L) polarized photons [101] 
\HM^)\ 2 = ^E e Hk 2 + ( 1 -^) 2 ] ejKKeflrM+mjKfcfm} (5.4) 

|^(^,r1,Q 2 )| 2 = l^^ e 2{ 4 Q2 z 2 (1 _^ ) 2^2 (e/ | fi | ) } ) (55) 

/ 

where a is the fine-structure constant, ej is the electric charge of the quark with flavor 
/, and Kq and K\ are the modified Bessel functions (McDonald functions). In the above 
expressions, 

e 2 = Zl {l- Zi )Q 2 + m) (5.6) 

controlls the transverse size(-distribution) of the emerging dipole, \fi\ oc 1/e/, that 
depends on the quark flavor through the current quark mass rrif. 

For small Q 2 , the perturbatively derived wave functions, (5.4) and (5.5), are not 
appropriate since the resulting color dipoles of size |r*j| oc 1/rrif S> lfm should en- 
counter non-perturbative effects such as confinement and chiral symmetry breaking. 
To take these effects into account the vector meson dominance (VMD) model [109] is 
usually used. However, the transition from the "partonic" behavior at large Q 2 to the 
"hadronic" one at small Q 2 can be modeled as well by introducing Q 2 -dependent quark 
masses, mf = mf(Q 2 ), that interpolate between the current quarks at large Q 2 and the 
constituent quarks at small Q 2 [52]. Following this approach, we use (5.4) and (5.5) 
also in the low-Q 2 region but with the running quark masses 

muAQ 2 ) = 0.178 GeV (1 - -^-) ®{Q 2 u4 - Q 2 ) , (5.7) 

o 2 

m s {Q 2 ) = 0.121 GeV + 0.129 GeV (1 - f^) B(Q 2 - Q 2 ) , (5.8) 
and the fixed charm quark mass 

m c = 1.25 GeV , (5.9) 

where the parameters Q 2 ld = 1-05 GeV 2 and Q 2 = 1.6 GeV 2 are taken directly from [52] 
while we reduced the values for the constituent quark masses rrif(Q 2 = 0) of [52] by 
about 20%. The smaller constituent quark masses are necessary to reproduce the total 
cross sections for j*p and 7*7* reactions at low Q 2 . Similar running quark masses are 
obtained in a QCD-motivated model of spontaneous chiral symmetry breaking in the 
instanton vacuum [110] that improve the description of 7*p scattering at low Q 2 [111]. 



5.2 Universal Energy Dependence 

Until now the T-matrix element (4.25) leads to energy independent total cross sections 
in contradiction to the experimental observation. In this section we introduce the energy 
dependence in a phenomenological way inspired by other successful models. 
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Most models for high-energy scattering are constructed to describe either hadron- 
hadron or photon-hadron reactions. For example, Kopeliovich et al. [112] as well as 
Berger and Nachtmann [57] focus on hadron-hadron scattering. In contrast, Golec- 
Biernat and Wusthoff [113] and Forshaw, Kerley, and Shaw [114] concentrate on photon- 
proton reactions. A model that describes the energy dependence in both hadron-hadron 
and photon-hadron reactions up to large photon virtualities is the two-pomeron model 
of Donnachie and Landshoff [62] . Based on Regge theory [35] , they find a soft pomeron 
trajectory with intercept 1 + e so f t ~ 1.08 that governs the weak energy dependence of 
hadron-hadron or j*P reactions with low Q 2 [115] and a hard pomeron trajectory with 
intercept 1 + ehard ~ 1-4 that governs the strong energy dependence of 7*p reactions 
with high Q 2 [62]. Similarly, we aim at a simultaneous description of hadron-hadron, 
photon-proton, and photon-photon reactions involving real and virtual photons as well. 

In line with other two-component (soft + hard) models [55,56,60,62,114] and the 
different hadronization mechanisms in soft and hard collisions, our physical ansatz de- 
mands that the perturbative and non-perturbative contributions do not interfere. There- 
fore, we modify the cosine-summation in (4.25) allowing only even numbers of soft and 
hard correlations, (x^) {x P ) with n, m £ IN. Interference terms with odd num- 
bers of soft and hard correlations are subtracted by the replacement 

cos [ex] = cos [c (x NP + X P )] cos [cx NP ] cos [cx P ] , (5.10) 

where c = 1/3 or 2/3. This prescription leads to the following factorization of soft and 
hard physics in the T-matrix element, 

T{s,t) = 2is Jd 2 b ± e 1 ^ fd Zl d 2 n fdz 2 d 2 r 2 |Vi^i,n)| 2 \M^r 2 )\ 2 

1 - H^ F ) cos (H - MH Ah")} ■ (5 - u> 

In the limit of small x-functions, x NP ^ 1 an d X P ^ 1) one gets 

T(s,t) = lis Jd 2 b ± e^ Jd Zl d 2 n jdz 2 d 2 r 2 |V>i(*i, n)| 2 |V 2 (^,r 2 )| 2 

*\[{X NP Y ' + (X P ) 2 } (5-12) 

so that the T-matrix element becomes a sum of a perturbative and a non-perturbative 
component. As already discussed in Sec. 4.3, the perturbative component (x P ) 2 co- 
incides with the well-known perturbative two-gluon exchange [99, 100] and the non- 
perturbative component (x^) 2 represents the corresponding non-perturbative gluonic 
interaction on the "two-gluon-exchange" level [24]. 

As the two-component structure of (5.12) reminds of the two-pomeron model of 
Donnachie and Landshoff [62], we adopt the powerlike energy increase and ascribe a 
weak energy dependence to the non-perturbative component x NP an d a strong one to 
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the perturbative component \ P 



(X NP ) 2 - (x^OO) 2 := (X NP ) 2 (~^) 

p 

(x p ) 2 - (^)) 2: =(^) 2 (^) e ( 5 - 13 ) 

with the scaling factor sqRq. The powerlike energy dependence with the exponents 
~ < e p < 1 guarantees Regge type behavior at moderately high energies, where 
the small- x limit (5.12) is appropriate. In (5.13) the energy variable s is scaled by the 
factor f 2 r 2 that allows to rewrite the energy dependence in photon-hadron scattering 
in terms of the appropriate Bjorken scaling variable x 

where \ f\\ is the transverse extension of the qq dipole in the photon. A similar factor 
has been used before in the dipole model of Forshaw, Kerley, and Shaw [114] and also 
in the model of Donnachie and Dosch [60] in order to respect the scaling properties 
observed in the structure function of the proton. 1 In the dipole-proton cross section of 
Golec-Biernat and Wiisthoff [113], Bjorken x is used directly as energy variable which 
is important for the success of the model. In fact, also in our model, the r 2 f 2 factor 
improves the description of 7*p reactions at large Q 2 . 

The powerlike Regge type energy dependence introduced in (5.13) is, of course, 
not mandatory but allows successful fits and can also be derived in other theoretical 
frameworks: A powerlike energy dependence is found for hard photon-proton reactions 
from the BFKL equation [116] and for hadronic reactions by Kopeliovich et al. [112]. 
However, these approaches need unitarization since their powerlike energy dependence 
will ultimately violate S'-matrix unitarity at asymptotic energies. In our model we use 
the following T-matrix element for investigations in the remaining chapters 

T(s,t) = 2is JcPb±e^ fd Zl d 2 n jdz 2 d 2 r 2 |Vi(«i,ri)| 2 |V>2(^,r 2 )| 2 (5.15) 
2 / 1 np / -A / 1 p / \\ 1 i 2 jvp /A / 2 p 



1 - 3 cos |^-x (s)j cos ^-x (s)J - - cos (s)J cos ^-x (s 

where the cosine functions ensure the unitarity condition in impact parameter space as 
shown in Chap. 6. Indeed, the multiple gluonic interactions associated with the higher 
order terms in the expansion of the cosine functions are important for the saturation 
effects observed within our model at ultra-high energies. 



1 In the model of Donnachie and Dosch [60], s\fi\ \n\ is used as the energy variable if both dipoles are 
small, which is in accordance with the choice of the typical BFKL energy scale but leads to discontinuities 
in the dipole-dipole cross section. In order to avoid such discontinuities, we use the energy variable (5.13) 
also for the scattering of two small dipoles. 
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Having ascribed the energy dependence to the ^-function, the energy behavior of 
hadron-hadron, photon-hadron, and photon-photon scattering results exclusively from 
the universal dipole-dipole scattering kernel. 

5.3 Model Parameters for High-Energy Scattering 

Lattice QCD simulations provide important information and constraints on the model 
parameters as discussed in Chap. 2. The fine tuning of the parameters was, however, 
directly performed on high-energy scattering data for hadron-hadron, photon-hadron, 
and photon-photon reactions [23] where an error (x 2 ) minimization was not feasible 
because of the non-trivial multi-dimensional integrals in the T-matrix element (5.15). 

The parameters a, k, G2, mo, M 2 , sqRq, e NP and e p determine the dipole-dipole 
scattering amplitude and are universal for all reactions considered. In addition, there 
are reaction-dependent parameters in the wave functions which are given in Sec. 5.1. 

The non-perturbative component involves the correlation length a, the gluon conden- 
sate G2, and the parameter k indicating the non-Abelian character of the correlator. 
With the simple exponential correlation function (2.50), we obtain the values given 
in (2.51) that have already been applied in Chap. 3 

a = 0.302 fm, k = 0.74, G 2 = 0.074 GeV 4 . 

The perturbative component involves the gluon mass mc as IR regulator (or inverse 
"perturbative correlation length" ) and the parameter M 2 that freezes the running cou- 
pling (2.47) for large distance scales. The following values are used in the following 
investigations 

ma = nip = 0.77 GeV and M 2 = 1.04 GeV 2 . (5.16) 

The energy dependence of the model is associated with the energy exponents e NP 
and e p , and the scaling parameter sqRq 

e NP = 0.125, e p = 0.73, and s R% = ( 47 GeV fm 2 ) 2 . (5.17) 

In comparison to the energy exponents of Donnachie and Landshoff [62, 115], e so ft ~ 0.08 
and ehard ~ 0.4, our exponents are significantly larger. However, the cosine functions in 
our T-matrix element (5.15) reduce the large exponents so that the energy dependence 
of the cross sections agrees with the experimental data as illustrated in Chap. 7. 
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Chapter 6 



Impact Parameter Profiles and 
Gluon Saturation 

In this chapter S-matrix unitarity constraints are considered in our model. On the basis 
of the impact parameter dependence of the scattering amplitude, saturation effects can 
be exposed that manifest the S'-matrix unitarity. For each impact parameter the energy 
at which the unitarity limit becomes important can be determined. The results are used 
to discuss gluon saturation and to localize saturation effects in experimental observables. 

The impact parameter dependence of the scattering amplitude is given by T(s, \ b±\), 

T(s,t = -q 2 ) = As jd 2 b ± e i ^T{s,\b 1 _\) , (6.1) 

and in particular by the profile function 

J(s,\b ± \) =2ImT( S ,|6J) , (6.2) 

which describes the blackness or opacity of the interacting particles as a function of 
the impact parameter |6j_| and the cm. energy ^/s. In fact, the profile function (6.2) 
determines all observables if the T-matrix is - as in our model - purely imaginary. 

6.1 S-Matrix Unitarity Constraints 

The S'-matrix unitarity, SS^ = S^S = H, leads directly to the unitarity condition in 
impact parameter space 1 [117, 118] 

JmT(s,\b ± \) = \T(s, \b ± \)\ 2 + G inel (s, \b ± \) , (6.3) 

integrating (6.3) over the impact parameter space and multiplying by a factor of 4 one obtains the 
relation a tot (s) = a el (s) + a mel (s). 



55 



where Gi ne i(s, \b±\) > is the inelastic overlap function [119]. This unitarity condition 
imposes an absolute limit on the profile function 

< 2|T( S ,|6J)| 2 < J(s,\b ± \) < 2 (6.4) 

and the inelastic overlap function, Gi ne i(s, \b±\) < 1/4. At high energies, the elastic 
amplitude is expected to be purely imaginary. Consequently, the solution of (6.3) reads 



J(S, \b ± \) = 1 ± V 1 - 4G -el( S ' ( 6 - 5 ) 

and leads with the minus sign corresponding to the physical situation to the reduced 
unitarity bound 

< J(s,|6l|) < 1 • (6.6) 

Reaching the black disc limit or maximum opacity at a certain impact parameter | b± | , 
■J( s j \b±\) = 1, corresponds to maximal inelastic absorption Gj ne ;(s, |6j_|) = 1/4 and equal 
elastic and inelastic contributions to the total cross section at that impact parameter. 

In our model every reaction is reduced to dipole-dipole scattering with well defined 
dipole sizes \fi\ and longitudinal quark momentum fractions Z{ as discussed in Sec. 4.1. 
Thus, the most general test of our model with respect to the unitarity constraints is 
performed with the profile function 

Jdd(s, \b±\,zi, |n |, z 2 , \r 2 \) = J Tj^jr 1 _ S DD (s,b±,z 1 ,r 1 ,z 2 ,f 2 ) , (6.7) 

where fa describes the dipole orientation, i.e. the angle between r*j and b±, and Sdd 
describes elastic dipole-dipole scattering 

" ' 4 NP t \ l f P / \\ , ^ ( NP i \ | f P/ 



Sdd = -co S ^-x^(s)) cos^-^( S )j + - cos ^- x ^(s)j cos ^( s )j (6.8) 

with the purely real-valued eikonal functions X NP ( S ) an d X P ( S ) defined in (5.13). Be- 
cause of |<Sd£>| < 1, a consequence of the cosine functions in (6.8) describing multiple 
gluonic interactions, Jdd respects the absolute limit (6.4). Thus, the elastic dipole- 
dipole scattering respects the unitarity condition (6.3). At high energies, the arguments 
of the cosine functions in Sdd become so large that these cosines average to zero in 
the integration over the dipole orientations. This leads to the black disc limit Jj^ff = 1 
reached at high energies first for small impact parameters. 



6.2 The Profile Function for Proton-Proton Scattering 

The profile function for proton-proton scattering 

J PP (s,\b±\) = Jdz 1 d 2 r 1 Jdz 2 d 2 r 2 \i{jp(z 1 , f 1 )\ 2 \ipp(z2, r 2 )\ 2 1 - S DD (s, b±, zi, fi, z 2 , f 2 

(6.9) 
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v/i = 10 1 GeV 
</s = 10 2 GeV 
y/s = 10 3 GeV 
y/s = 10 4 GeV 
y/s = 10 s GeV 
yfs~ = 10 6 GeV 
y/s~ = 10 7 GeV 
^ = 10 8 GeV 



\b ± \ [fm] 



Figure 6.1: The profile function for proton-proton scattering J pp (s, \b±\) as a function of the 
impact parameter |6j_| for cm. energies from y/s = 10 GeV to 10 s GeV. The unitarity limit (6.4) 
corresponds to J pp (s, \b±\) = 2 and the black disc limit (6.6) to J pp (s, \b±\) = 1. 



is obtained from (6.7) by weighting the dipole sizes |fj| and longitudinal quark momen- 
tum fractions Z{ with the proton wave function \ip p (zi, r*j)| 2 from Sec. 5.1. 

Using the model parameters from Sec. 5.3, we obtain the profile function J pp (s, \b±\) 
shown in Fig. 6.1 for cm. energies from y/s = 10 GeV to 10 8 GeV. Up to y/s ~ 100 GeV, 
the profile has approximately a Gaussian shape. Above y/s = 1 TeV, it develops into a 
broader and higher profile until the black disc limit is reached for y/s ~ 10 6 GeV and 
\b±\ = 0. At this point, the cosine functions in Sdd average to zero 

Jdz 1 d 2 r 1 Jdz 2 d 2 r 2 \ij p (z l ,r l )\ 2 \^ p (z2^ 2 )\ 2 S DD (yfs > 10 6 GeV, |6l| = 0, . . . ) w 

(6.10) 

so that the proton wave function normalization determines the maximum opacity 

J™- = j dzid 2 n j dz2( ? r2 |^(^,ri)| 2 \i; p (z 2 ,r 2 )\ 2 = 1 . (6.11) 

Once the maximum opacity is reached at a certain impact parameter, the height of the 
profile function saturates at that |&j_| while the width of the profile function extends 
towards larger impact parameters with increasing energy. Thus, multiple gluonic inter- 
actions important to respect the S'-matrix unitarity constraint (6.3) lead to saturation 
for y/s > 10 6 GeV. 
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The above behavior of the profile function illustrates the evolution of the proton 
with increasing cm. energy. The proton is gray and of small transverse size at small 
y/s but becomes blacker and more transversally extended with increasing y/s until it 
reaches the black disc limit in its center at y/s k, 10 6 GeV. Beyond this energy, the 
proton cannot become blacker in its central region but in its periphery with continuing 
transverse growth. Furthermore, the proton boundary stays diffusive as claimed also 
in [120]. 

According to our model the black disc limit will not be reached at LHC. Our predic- 
tion of y/s f» 10 6 GeV = 10 3 TeV for the onset of the black disc limit in proton-proton 
collisions is about two orders of magnitude beyond the LHC energy y/s = 14 TeV. This 
is in contrast to, e.g. [121], where the value predicted for the onset of the black disc limit 
is y/s = 2 TeV, i.e. small enough to be reached at LHC. However, note that our profile 
function J pp (s, \ b±\) yields good agreement with experimental data for cross sections up 
to the highest energies as shown in Chap. 7. 

For hadron-hadron reactions in general, the wave function normalization of the 
hadrons determines the maximum opacity analogous to (6.11) and the transverse hadron 
size the cm. energy at which it is reached. Consequently, the maximum opacity obtained 
for np and Kp scattering is identical to the one for pp scattering due to the normal- 
ization (5.2). The smaller size of pions and kaons in comparison to protons, however, 
demands slightly higher cm. energies to reach this maximum opacity. Such size effects 
become more pronounced in longitudinal photon-proton scattering, where the size of 
the dipole emerging from the photon can be controlled by the photon virtuality. 



6.3 The Profile Function for Photon-Proton Scattering 

The profile function for a longitudinal photon scattering off a proton p 

J^ L p{s,\b±\,Q 2 ) = Jdz 1 d 2 r 1 Jdz 2 d 2 r 2 \iJ 1 i{zi,r 1 ,Q 2 )\ 2 \i) p (z2,r2)\ 2 

1 - S DD (bj_,s,z 1 ,f 1 ,Z2,f 2 ) (6.12) 

is calculated with the longitudinal photon wave function |^ 7 £ (zi, fi, Q 2 )\ 2 given in (5.5). 
In this way, the profile function (6.12) is ideally suited for the investigation of dipole 
size effects since the photon virtuality Q 2 determines the transverse size of the dipole 
into which the photon fluctuates before it interacts with the proton. 

Figure 6.2 shows the dependence of the profile function J 7 £ p (s, \b±\, Q 2 ) divided 
by a/ir at a photon virtuality of Q 2 = 1 GeV 2 for cm. energies y/s from 10 GeV to 
10 9 GeV, where a is the fine-structure constant. One clearly sees that the qualitative 
behavior of this rescaled profile function is similar to the one for proton-proton scatter- 
ing. However, the black disc limit induced by the underlying dipole-dipole scattering 
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■s/a = 10 1 GeV 
yCs = io 2 GeV 
v/i = 10 3 GeV 
y/a = 10 4 GeV 
= 10 s GeV 
^1 = io 6 GeV 
v^i = IO 7 GeV 
V 7 ! = 10 8 GeV 
Vs = 10 9 GeV 



[fill] 



Figure 6.2: The profile function for a longitudinal photon scattering off a proton J 7 ^ p (s, \b±\, Q 2 ) 
divided by ol/tt as a function of the impact parameter at a photon virtuality of Q 2 = 1 GeV 2 
and cm. energies from y/s = 10 GeV to 10 9 GeV. The value of the black disc limit is J™l a *(Q 2 = 



lGeV 



0.00164 



depends on the photon virtuality Q 2 and is given by the normalization of the longitu- 
dinal photon wave function 

J^(Q 2 ) = Jdzd 2 r\i; rL (z,f,Q 2 )\ 2 (6.13) 

since the proton wave function is normalized to one. 

The photon virtuality Q 2 does not only determine the absolute value of the black 
disc limit but also the cm. energy at which it is reached. This is illustrated in Fig. 6.3, 
where the yfs dependence of J 7 £ p (s, |6j_| = 0, Q 2 ) divided by a/ir is presented for Q 2 = 
1, 10, and 100 GeV 2 . With increasing resolution Q 2 , i.e. decreasing dipole sizes, I^VE | 2 oc 
1/Q 2 , the absolute value of the black disc limit grows and higher energies are needed 
to reach this limit. 2 The growth of the absolute value of the black disc limit is simply 
due to the normalization of the longitudinal photon wave function while the need for 
higher energies to reach this limit is due to the decreasing interaction strength with 
decreasing dipole size. The latter explains also why the energies needed to reach the 
black disc limit in irp and Kp scattering are higher than in pp scattering. Comparing 

2 Note that the Bjorken x at which the black disc limit is reached decreases with increasing photon 
virtuality Q 2 . See also Fig. 6.4. 
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Figure 6.3: The profile function for a longitudinal photon scattering off a proton J 7 * p (s, Q 2 ) 
divided by a/n as a function of the cm. energy y/s at zero impact parameter = 0) for 
photon virtualities of Q 2 = 1, 10, and 100 GeV 2 . 



j^P scattering at Q 2 = 1 GeV 2 with pp scattering quantitatively, the black disc limit 
jmax(g2 _ 1Ge V 2 ) = 0.00164 is about three orders of magnitude smaller because of 

the photon wave function normalization (oc a /it). At \b±\ = it is reached at an energy 
of y/s ~ 10 s GeV, which is about two orders of magnitude higher because of the smaller 
dipoles involved. 

The way in which the profile function J 7 £ p (s, \b±\, Q 2 ) approaches the black disc 
limit at high energies depends on the shape of the proton and longitudinal photon wave 
function at small dipole sizes |ri a|. At high energies, dipoles of typical sizes < |ri 2I < 
i?o (so/s) 1 ^ give the main contribution to S 7 * p = 1 — J 7 £ P because of (5.13) and the 
fact that the contribution of the large dipole sizes averages to zero upon integration 
over the dipole orientations, cf. also (6.10). Since 5* 7 £ P is a measure of the proton 
transmittance, this means that only small dipoles can penetrate the proton at high 
energies. Increasing the energy further, even these small dipoles are absorbed and the 
black disc limit is reached. However, the dependence of the profile function on the short 
distance behavior of normalizable wave functions is weak which can be understood as 
follows. Because of color transparency, the eikonal functions X ( s ) an d X P ( S ) are small 
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for small dipole sizes < {fi^l < Ro (so/s) 1 ^ at large s/s. Consequently, Sdd ~ 1 and 

1 r c (s) 1 r c (s) 

« J^(Q 2 )-4ir 2 Jd Zl Jdnnl^iz^r^Q 2 )^ fdz 2 j dr 2 r 2 \^ p (z 2 ,r 2 )\ 2 (6.14) 


where r c (s) ~ i?o (so/s) 1 ^ 4 . Clearly, the linear behavior from the phase space factors r± t2 
dominates over the ri^-dependence of normalizable wave functions. 3 More generally, 
for any profile function involving normalizable wave functions, the way in which the 
black disc limit is approached depends only weakly on the short distance behavior of 
the wave functions. 



6.4 A Scenario for Gluon Saturation 



In this section we discuss saturation of the impact parameter dependent gluon distri- 
bution of the proton xG(x,Q 2 , \b±\). Using a leading twist, next-to-leading order QCD 
relation between xG(x, Q 2 ) and the longitudinal structure function Fl(x, Q 2 ), we relate 
xG(x, Q 2 , \b±\) to the profile function J 7 £ p (s = Q 2 /x, \b±\, Q 2 ) and find low-x saturation 

of xG(x, Q 2 , \ b±\) as a manifestation of S'-matrix unitarity. The resulting xG(x, Q 2 , \ b±\) 
is, of course, only an estimate since our profile function contains also higher twist contri- 
butions. Furthermore, in the considered low-x region, the leading twist, next-to-leading 
order QCD formula may be inadequate as higher twist contributions [122] and higher 
order QCD corrections [123, 124] are expected to become important. Nevertheless, still 
assuming a close relation between Fl(x, Q 2 ) and xG(x, Q 2 ) at low x, we think that our 
approach provides some insight into the gluon distribution as a function of the impact 
parameter and into its saturation. 

The gluon distribution of the proton xG(x,Q 2 ) is defined as follows: xG(x,Q 2 )dx 
gives the momentum fraction of the proton which is carried by the gluons in the in- 
terval [x,x + dx] as seen by probes of virtuality Q 2 . The impact parameter dependent 
gluon distribution xG(x, Q 2 , \b±\) is the gluon distribution xG(x, Q 2 ) at a given impact 
parameter |6j_| so that 

xG(x,Q 2 ) = Jd 2 b ± xG(x,Q 2 ,\b ± \) . (6.15) 

In leading twist, next-to-leading order QCD, the gluon distribution xG(x,Q 2 ) is 

3 For our choice of the wave functions in (6.14), one sees very explicitly that the specific Gaussian 
behavior of \tp p (z2, r*2 ) | 2 and the logarithmic short distance behavior of \tp~/i (zi, r\, Q 2 )\ 2 is dominated 
by the phase space factors ri j2 - 
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related to the structure functions Fl(x, Q 2 ) and ^(x, Q 2 ) of the proton [125] 



Fl(x,Q 2 ) 



a. 



7T 



x y \yj V Jx v 



l dy( 



x 

\y 



yG(y,Q 2 ) 



(6.16) 

where J2f e / i s a flavor sum over the quark charges squared. For four active flavors and 
x < 10~ 3 , (6.16) can be approximated as follows [126] 



xG(x,Q 2 ) -5.8 



37T 1 

— F L (0Al7x, Q 2 ) - — F 2 (0.75x, Q 2 ) 



(6.17) 



For typical A QC d = 100 - 300 MeV and Q 2 = 50 - 100 GeV 2 , the coefficient of F L in 
(6.17), 37r/(4a s ) = O(10), is large compared to the one of Fi- Taking into account 
also the values of F2 and Fl, the gluon distribution is mainly determined by the lon- 
gitudinal structure function for x < 10~ 3 in this Q 2 region. The longitudinal structure 
function can be expressed in terms of the profile function for longitudinal photon-proton 
scattering using the optical theorem (cf. (7.1)) 



F L (x,Q 2 



Q 2 



47r 2 a 



Q 2 



in 2 a 



2 ld 2 b ± J^ LP (x,\b ± \,Q 2 



(6.18) 



where the s-dependence of the profile function is rewritten in terms of the Bjorken 
scaling variable x = Q 2 /s. Thus, neglecting the F2 term in (6.17), the gluon distribution 
reduces to 

xG(x,Q 2 ) « 1.305 ■ U2u 1 i»- 1 ^ 



7r 2 a s a 



- d z b ± J llp (0A17x,\b ± \,Q 2 



(6.19) 



Comparing (6.15) with (6.19), it seems natural to relate the integrand of (6.19) to the 
impact parameter dependent gluon distribution 



O 2 TV 

xG(x,Q 2 ,\b ± \) « 1.305 -f J 7£p (0.417x,|6 ± |,Q 2 ) . 

7T GLg (y. 



(6.20) 



The black disc limit of the profile function for longitudinal photon-proton scatter- 
ing (6.13) imposes accordingly an upper bound on xG(x,Q 2 , \b_\_\) 

xG(x, Q 2 , \b ± \) < xG max (Q 2 ) « 1.305 - J^{Q 2 ) , (6.21) 

which is the low-x saturation value of the gluon distribution xG(x, Q 2 , \b±\) in our 
approach. With vrJ™p X (Q 2 )/a « 1 as shown in Fig. 6.3, a compact approximation 
of (6.21) is obtained 

xG(x,Q 2 ,\b ± \) < xG max (Q 2 ) « 4- > ( 6 - 22 ) 
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Figure 6.4: The gluon distribution of the proton at zero impact parameter xG(x, Q 2 , \b±\ = 0) 
as a function of x for Q 2 — 1, 10, andlOOGeV 2 . The results are obtained within approxima- 
tion (6.20). 



which is consistent with the results in [124, 127, 128] and indicates strong color field 
strengths Cr" ~ \j ^fa~ s as well. 

According to our relations (6.20) and (6.21), the blackness described by the profile 
function is a measure for the gluon distribution and the black disc limit corresponds to 
the maximum gluon distribution reached at the impact parameter under consideration. 
In accordance with the behavior of the profile function J 7 £ P , see Fig. 6.2, the gluon 
distribution xG(x, Q 2 , \ b±\) decreases with increasing impact parameter for given values 
of x and Q 2 . Consequently, the gluon density has its maximum in the geometrical center 
of the proton, i.e. at zero impact parameter, and decreases towards the periphery. With 
decreasing x at given Q 2 , the gluon distribution xG(x, Q 2 , \b±\) increases and extends 
towards larger impact parameters just as the profile function J^* p for increasing s. The 
saturation of the gluon distribution xG{x, Q 2 , \ b±\) sets in first in the center of the proton 
(|6j_| =0) at very small Bjorken x. 

In Fig. 6.4 the gluon distribution xG(x, Q 2 , \b±\ = 0) is shown as a function of x 
for Q 2 = 1, 10, andlOOGeV 2 , where relation (6.20) has been used also for low photon 
virtualities. Evidently, the gluon distribution xG(x, Q 2 , \b±\ = 0) saturates at very 
low values of x < 10~ 10 for Q 2 > lGeV 2 . The photon virtuality Q 2 determines the 
saturation value (6.21) and the Bjorken x at which it is reached; cf. also Fig. 6.2. For 
larger Q 2 , the low-x saturation value is larger and is reached at smaller values of x, as 
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claimed also in [129]. Moreover, the growth of xG(x, Q 2 , \b±\ = 0) with decreasing x 
becomes stronger with increasing Q 2 . This results from the stronger energy increase of 
the perturbative component, e p = 0.73, that becomes more important with decreasing 
dipole size. 

According to our approach, the onset of the xG(x,Q 2 , |&i|)-saturation appears for 
Q 2 > lGeV 2 at x < KT 10 , which is far below the x-region accessible at HERA (x > 
10~ 6 ). Even for THERA (x > 10" 7 ), gluon saturation is not predicted for Q 2 > 1 GeV 2 . 
However, since the HERA data can be described by models with and without saturation 
embedded [129], the present situation is not conclusive. 4 

Note that the S'-matrix unitarity condition (6.3) together with (6.20) requires the 
saturation of the impact parameter dependent gluon distribution xG(x, Q 2 , \b±\) but not 
the saturation of the integrated gluon distribution xG(x,Q 2 ). Due to multiple gluonic 
interactions in our model, this requirement is fulfilled, as can be seen from Fig. 6.2 
and relation (6.20). Indeed, approximating the gluon distribution xG(x, Q 2 , \ b±\) in the 
saturation regime of very low x by a step-function 

xG(x,Q 2 ,\b±\) « xG max (Q 2 ) Q(R(x,Q 2 ) - \b±\) , (6.23) 

where R(x,Q 2 ) denotes the full width at half maximum of the profile function, one 
obtains with (6.15), (6.21) and (6.22) the integrated gluon distribution 

xGix , Q3) M L305 1 J^tft) » Q2R2(X ' Q2) . (0.24) 

na s a ira s 

which does not saturate because of the increase of the effective proton radius R(x, Q 2 ) 
with decreasing x. Nevertheless, although xG(x,Q 2 ) does not saturate, the saturation 
of xG(x, Q 2 , \b±\) leads to a slow-down in its growth towards small x. 5 Interestingly, 
our result (6.24) coincides with the result of Mueller and Qiu [124]. 

Finally, it must be emphasized that the low-x saturation of xG(x, Q 2 , \ b±\), required 
in our approach by the S'-matrix unitarity, is realized by multiple gluonic interactions. 
In other approaches that describe the evolution of the gluon distribution with varying 
x and Q 2 , gluon recombination leads to gluon saturation [123,124,131-133], which is 
reached when the probability of a gluon splitting up into two is equal to the probability 
of two gluons fusing into one. A more phenomenological understanding of saturation is 
attempted in [113, 134]. 



4 So far, the most striking hint for saturation in the present HERA data at x ~ 1CP 4 and Q 2 < 2 GeV 2 
has been the turnover of dFiix, Q 2 )/dln(Q 2 ) towards small x in the Caldwell plot [130], which is still 
a controversial issue due to the correlation of Q 2 and x values. 

This is analogous to the rise of the total pp cross section with growing cm. energy that slows down 
as the corresponding profile function J pp (s, \b±\) reaches its black disc limit as shown in Sec. 7.1. 
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Chapter 7 



Comparison with Experimental 
Data 



In this chapter we present the phenomenological performance of our model. We compute 
total, differential, and elastic cross sections, structure functions, and diffractive slopes 
for hadron-hadron, photon-proton, and photon-photon scattering, compare the results 
with experimental data including cosmic ray data, and provide predictions for future 
experiments. Having studied the saturation of the impact parameter profiles, we show 
here how this manifestation of unitarity translates into the quantities mentioned above 
and how it could become observable. 

Using the T-matrix (5.15) with the wave functions and parameters from Sees. 5.1 
and 5.3, we compute the pomeron contribution to pp, pp, 7r ± p, K^p, j*p, and 77 
reactions in terms of the universal dipole-dipole scattering amplitude Sdd- This allows 
one to compare reactions induced by hadrons and photons in a systematic way. In fact, 
it is our aim to provide a unified description of all these reactions and to show in this 
way that the pomeron contribution to the above reactions is universal and can be traced 
back to the dipole-dipole scattering amplitude Sdd- 

Our model describes pomeron (C = +1 gluon exchange) but neither odderon (C = 
— 1 gluon exchange) nor reggeon exchange (quark-antiquark exchange) as discussed in 
Sec. 4.3. Only in the computation of the hadronic total cross sections the reggeon 
contribution is added [59, 115]. This improves the agreement with the data for y/s < 
100 GeV and describes exactly the differences between ab and ah reactions. 

The fine tuning of the model and wave function parameters was performed on the 
data shown below. The resulting parameter set given in Sees. 5.1 and 5.3 is used 
throughout this chapter. 



65 



7.1 Total Cross Sections 



The total cross section for the high-energy reaction ab — > X is related via the optical 
theorem to the imaginary part of the forward elastic scattering amplitude and can also 
be expressed in terms of the profile function (6.2) 

o%(s) = - s lmT(s,t = 0) = 2 Jd 2 b ± J ab (s,\b ± \) , (7.1) 

where a and b label the initial particles whose masses were neglected as they are small 
in comparison to the cm. energy yfs. 

We compute the pomeron contribution to the total cross section, cr 1 "^ (s), from the 
T-matrix (5.15), as explained above, and add only here a reggeon contribution of the 
form [59, 115] 

/ \ -0.4525 

= (j^f) > (7-2) 

where X ab depends on the reaction considered: X pp = 56.08 mb, X p p = 98.39 mb, 
X w + p = 27.56 mb, X w - p = 36.02 mb, X K + p = 8.15 mb, X K - p = 26.36 mb, X Jp = 
0.129 mb, and X 77 = 605 nb. Accordingly, we obtain the total cross section 

<>%{*) = +*%"{s) (7.3) 

for pp, pp, ir^p, K^p, 7p, and 77 scattering. 

The good agreement of the computed total cross sections with the experimental 
data is shown in Fig. 7.1. Here the solid lines represent the theoretical results for pp, 
ir + p, K + p, jp, and 77 scattering and the dashed lines the ones for pp, n~p, and K~p 
scattering. The pp, pp, ir^p, K^p, ^p [10] and 77 data [135] taken at accelerators 
are indicated by the closed circles while the closed squares (Fly's eye data) [136] and 
the open circles (Akeno data) [137] indicate cosmic ray data. Only real photons are 
considered which are, of course, transverse polarized. 

The prediction for the total pp cross section at LHC (yfs = 14TeV) is (7*°* = 
114.2 mb in good agreement with cosmic ray data. Compared with other works, our 
LHC prediction is close to the one of Block et al. [138], a p ° p = 108 ± 3.4 mb, but 
considerably larger than the one of Donnachie and Landshoff [115], cr*°* = 101.5 mb. 

The differences between ab and ab reactions for y/s < 100 GeV result solely from 
the different reggeon contributions which die out rapidly as the energy increases. The 
pomeron contribution to ab and ab reactions is, in contrast, identical and increases as 
the energy increases. It thus governs the total cross sections for y/s > 100 GeV where 
the results for ab and ab reactions coincide. 

The differences between pp (pp), ^p, and K^p scattering result from the different 
transverse extension parameters, S p = 0.86 fm > S n = 0.607 fm > Sk = 0.55 fm, 
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Figure 7.1: The total cross section a tot as a function of the cm. energy yfs for pp, pp, ir^p, 
K p, 7p, and 77 scattering. The solid lines represent the model results for pp, ir + p, K + p, jp 
and 77 scattering and the dashed lines the ones for pp, ir~p, and K~p scattering. The pp, pp, 
n^p, K ± p, ryp [10] and 77 data [135] taken at accelerators are indicated by the closed circles 
while the closed squares (Fly's eye data) [136] and the open circles (Akeno data) [137] indicate 
cosmic ray data. The arrows at the top point to the LHC energy, s/s = 14 TeV, and to the onset 
of the black disc limit in pp (pp) reactions, * 10 6 GeV. 
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cf. Sec. 5.1. Since a smaller transverse extension parameter favors smaller dipoles, 
the total cross section becomes smaller and the short distance physics described by 
the perturbative component becomes more important which leads to a stronger energy 
growth due to e p = 0.73 > e NP = 0.125. In fact, the ratios a** /a™ and a^/a^ 
converge slowly towards unity with increasing energy as can already be seen in Fig. 7.1. 

For real photons, the transverse size is governed by the constituent quark masses 
nif(Q 2 = 0), cf. Sec. 5.1, where the light quarks have the strongest effect, i.e. at^ oc 
and oc l/m^j. Furthermore, in comparison with hadron-hadron scatter- 
ing, there is the additional suppression factor of a for 7p and a 2 for 77 scattering 
coming from the photon-dipole transition. In the 77 reaction, also the box diagram 
contributes [59, 139] but is neglected since its contribution to the total cross section is 
less than 1% [60]. 

It is worthwhile mentioning that total cross sections for pp (pp) , 7r ± p, and K ± p scat- 
tering do not depend on the longitudinal quark momentum distribution in the hadrons 
since the underlying dipole-dipole cross section is independent of the longitudinal quark 
momentum fraction Zi for t = 0. We have shown this analytically on the two-gluon- 
exchange level in [24]. 

Saturation effects as a manifestation of the S'-matrix unitarity can be seen in Fig. 7.1. 
Having investigated the profile function for pp (pp) scattering, we know that this profile 
function becomes higher and broader with increasing energy until it saturates the black 
disc limit first for zero impact parameter = 0) at y/s ~ 10 6 GeV. Beyond this 
energy, the profile function cannot become higher but expands towards larger values of 
|&i_|. Consequently, the total cross section (7.1) increases no longer due to the growing 
blackness at the center but only due to the transverse expansion of the hadrons. This 
tames the growth of the total cross section as can be seen for the total pp cross section 
beyond ^ ~ 10 6 GeV in Fig. 7.1. 

At energies beyond the onset of the black disc limit at zero impact parameter, the 
profile function can be approximated by 

JT OX (s, \b ± \) = N a N b e(R(s) - \b ± \) (7.4) 

where N a ^ denotes the normalization of the wave functions of the scattered particles 
and R(s) the full width at half maximum of the exact profile function J a b(s, \ b±_\) that 
reflects the effective radii of the interacting particles. Thus, the energy dependence 
of the total cross section (7.1) is driven exclusively by the increase of the transverse 
extension of the particles R(s) 

a^(s) = 2irN a N b R( S ) 2 , (7.5) 

which is known as geometrical scaling [118, 140]. The growth of R(s) can at most be 
logarithmic for \/s — > 00 because of the Froissart bound [18]. In fact, a transition from 
a power-like to an ln 2 -increase of the total pp cross section seems to set in at about 
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yfs ~ 10 6 GeV as visible in Fig. 7.1. Moreover, since the hadronic cross sections join 
for ^/s — ► oo, becomes independent of the hadron-hadron reaction considered at 
asymptotic energies as long as N a ^ = 1. Also for photons of different virtuality Q\ and 
Q\ one can check that the ratio of the total cross sections o^i^Q 2 ) j oi°i v (Q\) converges 
to unity at asymptotic energies in agreement with the conclusion in [141]. 



7.2 The Proton Structure Function 



The total cross section for the scattering of a transverse (T) and longitudinally (L) 
polarized photon off the proton, rf*, LP {%, Q 2 ), at photon virtuality Q 2 and cm. energy 1 

squared s = Q 2 /x is equivalent to the structure functions of the proton 

%(x,0 2 ) = ^ iP (x,Q 2 ) (7.6) 



and 

F 2 (x, Q 2 ) = F T {x, Q 2 ) + F L (x, Q 2 ) . (7.7) 

Reactions induced by virtual photons are particularly interesting because the trans- 
verse separation of the quark-antiquark pair that emerges from the virtual photon de- 
creases as the photon virtuality increases (cf. Sec. 5.1) 

(7.8) 



^Q 2 + 4<d ' 

where ra u ^ is a mass of the order of the constituent u-quark mass. With increasing 
virtuality, one probes therefore smaller transverse distance scales of the proton. 

In Fig. 7.2 the Q 2 -dependence of the total j*p cross section 

o*? p (s, Q 2 ) = o*g p (s, Q 2 ) + a% p (s, Q 2 ) (7.9) 

is presented, where the model results (solid lines) are compared with the experimental 
data for cm. energies from y/s = 20 GeV up to 245 GeV. Note the indicated scaling 
factors at different yfs values. The low-energy data at yfs = 20 GeV are from [142] 
while the data at higher energies have been measured at HERA by the HI [11, 143] and 
ZEUS collaboration [13, 144]. At Q 2 = 0.012 GeV 2 , also the photoproduction (Q 2 = 0) 
data from [145] are displayed. 

The model results are in reasonable agreement with the experimental data in the 
window shown in Fig. 7.2. The total 7*p cross section levels off towards small values 
of Q 2 as soon as the photon size |r 7 |, i.e the resolution scale, becomes comparable to 



Here yfs refers to the cm. energy in the 7*p system. 
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Figure 7.2: The total 7*p cross section <r*°p(s, Q 2 ) as a function of the photon virtuality Q 2 
for cm. energies from yfs — 20 GeV to 245 GeV, where the model results (solid lines) and the 
experimental data at different \/s values are scaled with the indicated factors. The low energy 
data at y/s = 20 GeV are from [142], the data at higher energies from the HI [11,143] and 
ZEUS collaboration [13,144]. The photoproduction (Q 2 = 0) data from [145] are displayed at 
Q 2 = 0.012 GeV 2 . 
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the proton size. Our model reproduces this behavior by using the perturbative photon 
wave functions with Q 2 -dependent quark masses, mf(Q 2 ), that interpolate between the 
current (large Q 2 ) and the constituent (small Q 2 ) quark masses as explained in detail 
in Sec. 5.1. The decrease of cr*°* p with increasing Q 2 results from decreasing dipole sizes 
and the fact that small dipoles do not interact as strongly as large dipoles. 

The x-dependence of the computed proton structure function Fz(x, Q 2 ) is shown in 
Fig. 7.3 for Q 2 = 0.3, 2.5, 12, and 120 GeV 2 in comparison to the data measured by the 
HI [146] and ZEUS [147] detector. Within our model, the increase of -^(x, Q 2 ) towards 
small Bjorken x becomes stronger with increasing Q 2 in agreement with the trend in 
the HERA data. This behavior results from the fast energy growth of the perturbative 
component that becomes more important with increasing Q 2 due to the smaller dipole 
sizes involved. 
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Figure 7.3: The x-dependence of the computed proton structure function i^x, Q 2 ) (solid line) 
for Q 2 = 0.3, 2.5, 12, and 120 GeV 2 in comparison to the data measured by the HI [146] and 
ZEUS [147] detector, and the low energy data at ^/s = 20 GeV from [142]. 

As can be seen in Fig. 7.3, the data show a faster increase with decreasing x than the 
model outside the low-Q 2 region. This results from the weak energy boost of the non- 
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perturbative component that dominates F2(x,Q 2 ) in our model. In fact, even for large 
Q 2 the non-perturbative contribution overwhelms the perturbative one, which explains 
also the overestimation of the data for x > 10^ 3 . 

This problem is typical for the SVM model applied to the scattering of a small size 
dipole off a proton. In an earlier application by Riiter [55], an additional cut-off was 
introduced to switch from the non-perturbative to the perturbative contribution as soon 
as one of the dipoles is smaller than r cut = 0.16 fm. This yields a better agreement with 
the data also for large Q 2 but leads to a discontinuous dipole-proton cross section. In 
the model of Donnachie and Dosch [60], a similar SVM-based component is used also for 
dipoles smaller than R c = 0.22 fm with a strong energy boost instead of a perturbative 
component. Furthermore, their SVM-based component is tamed for large Q 2 by an 
additional a s {Q 2 ) factor. 

We did not follow these lines in order to keep a continuous, Q 2 -independent dipole- 
proton cross section and, therefore, cannot improve the agreement with the F2(x,Q 2 ) 
data without losing quality in the description of hadronic observables. Since our non- 
perturbative component relies on lattice QCD, we are more confident in describing 
non-perturbative physics and thus put more emphasis on the hadronic observables. 
Admittedly, our model misses details of the proton structure that become visible with 
increasing Q 2 . However, most other existing models provide neither the profile functions 
nor a simultaneous description of hadronic and 7*-induced processes. 



7.3 The Slope B of Elastic Forward Scattering 



The local slope of elastic scattering B(s,t) is defined as 



B(M) : = S Un 



do- 



el 



alt 



■(s,t) 



(7.10) 



and characterizes the diffractive peak of the differential elastic cross section da el /dt(s, t) 
discussed below. Here we concentrate on the slope for elastic forward (t = 0) scattering 
also called slope parameter 



B(s) :=B(s,t = 0) = - 



lfd 2 b ± \b ± \ 2 J(s,\b ± \)_ = l {b2) 



2 Jd 2 b ± J(s,\b ± 



(7.11) 



which measures the root mean squared interaction radius (b 2 ) of the scattered particles, 
and does not depend on the opacity. 

We compute the slope parameter with the profile function from the T-matrix (5.15) 
and neglect the reggeon contributions, which are relevant only at small cm. energies, 
so that the same result is obtained for ab and ah scattering. 
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Figure 7.4: The elastic slope parameter B(s) as a function of the cm. energy y/s for pp and pp 
forward (t = 0) scattering. The solid line represents the model result that is compared to the 
data for pp (open circles) and pp (closed circles) reactions from [148-150]. 



In Fig. 7.4 the resulting slope parameter B(s) is shown as a function of yfs for pp and 
pp scattering (solid line) and compared with the pp (open circles) and pp (closed circles) 
data from [148-150]. As expected from the opacity independence of the slope parameter 
(7.11), saturation effects as seen in the total cross sections do not occur. Indeed, one 
observes an approximate B(s) oc R 2 (s) oc \r? {yfs / ^J~sq) growth for yfs > 10 4 GeV. This 
behavior agrees, of course, with the transverse expansion of J pp (s, \b±\) for increasing 
y/s shown in Fig. 6.1. Analogous results are obtained also for np and Kp scattering. 

For the good agreement of our model with the data, the finite width of the longitudi- 
nal quark momentum distribution in the hadrons, i.e. Az p , Az-^, and Azk ^ in (5.1), 
is important as the numerator in (7.11) depends on this width. In fact, B(s) comes out 
more than 10% smaller with Az p , Az n , and Azx = 0. Furthermore, a strong growth of 
the perturbative component, e p = 0.73, is important to achieve the increase of B(s) for 
yfs > 500 GeV indicated by the data. 

It must be emphasized that only the simultaneous fit of the total cross section and the 
slope parameter provides the correct shape of the profile function J(s, This shape 
leads then automatically to a good description of the differential elastic cross section 
da el /dt(s,t) as demonstrated below. Astonishingly, only few phenomenological models 
provide a satisfactory description of both quantities [112, 138]. In the approach of [57], 
for example, the total cross section is described correctly while the slope parameter 
exceeds the data by more than 20% already at yfs = 23.5 GeV and thus indicates 
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deficiencies in the form of J(s, \b±\). 



7.4 The Differential Elastic Cross Section 



The differential elastic cross section obtained from the squared absolute value of the 
T-matrix element 

da el 1 

= iis?i^'>P (7 - 12> 

can be expressed for our purely imaginary T-matrix (5.15) in terms of the profile function 

da el , , 1 



(7.13) 



and is, thus, very sensitive to the transverse extension and opacity of the scattered 
particles. Equation (7.13) shows the analogy to optical diffraction, where J(s, \b±\) 
describes the distribution of an absorber that causes the diffraction pattern observed 
for incident plane waves. 

In Fig. 7.5 the differential elastic cross section computed for pp and pp scattering 
(solid line) is shown as a function of \t\ = cf 2 at yfs = 23.5, 30.7, 44.7, 63, 546, and 
1800 GeV and compared with experimental data (open circles), where the pp data at 
y/s = 23.5, 30.7, 44.7, and 63 GeV were measured at the CERN ISR [140], the pp data 
at y/s = 546 GeV at the CERN SppS [149], and the pp data at y/s = 1.8 TeV at the 
Fermilab Tevatron [150, 151]. The prediction of our model for the pp differential elastic 
cross section at the CERN LHC, y/s = 14 TeV, is given in Fig. 7.6. 

For all energies, the model reproduces the experimentally observed diffraction pat- 
tern, i.e the characteristic diffraction peak at small \t\ and the dip structure at medium 
\t\. As the energy increases, also the shrinking of the diffraction peak is described 
which reflects the rise of the slope parameter B(s,t = 0) already discussed above. The 
shrinking of the diffraction peak comes along with a dip structure that moves towards 
smaller values of \t\ as the energy increases. This motion of the dip is also described 
approximately. 

The dip in the theoretical curves reflects a change of sign in the T-matrix ele- 
ment (5.15). As the latter is purely imaginary, it is not surprising that there are devia- 
tions from the data in the dip region. Here the real part is expected to be important [151] 
which is in the small |t| region negligible in comparison to the imaginary part. 

The difference between the pp and pp data, a deep dip for pp but only a bump 
or shoulder for pp reactions, requires a C = — 1 contribution. Besides the reggeon 
contribution at small energies, 2 one expects here an additional perturbative C = —I 

2 Zooming in on the result for y/s — 23.5 GeV, one finds an underestimation of the data for all \t\ 
before the dip, which is correct as it leaves room for the reggeon contribution being non-negligible at 
small energies. 
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Figure 7.5: The differential elastic cross section for pp and pp scattering as a function of \t\. The 
result of our model, indicated by the solid line, is compared for ^/s = 23.5, 30.7, 44.7, and 63 GeV 
with the CERN ISR pp data [140], for = 546 GeV with the CERN SppS data [149], and for 
= 1.8 TeV with the Fermilab Tevatron pp data [150, 151], all indicated by the open circles 
with error bars. 
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Figure 7.6: The prediction of our model for the pp differential elastic cross section at LHC 
(^/s = 14 TeV) as a function of momentum transfer \t\ up to 1 GeV 2 . 

contribution such as three-gluon exchange [152,153] or an odderon [154-156]. In fact, 
allowing a finite size diquark in the (anti-)proton an odderon appears that supports the 
dip in pp but leads to the shoulder in pp reactions [156]. 

For the differential elastic cross section at the LHC energy, ^/s = 14 TeV, the above 
findings suggest an accurate prediction in the small-|i| region but a dip at a position 
smaller than the predicted value at \t\ ~ 0.35 GeV 2 . Our confidence in the validity of 
the model at small |t| is supported additionally by the total cross section that fixes 
da el /dt(s,t = 0) and is in agreement with the cosmic ray data shown in Fig. 7.1. 
Concerning our prediction for the dip position, it is close to the value \t\ ~ 0.41 GeV 2 
of [138] but significantly below the value \t\ ~ 0.55 GeV 2 of [57]. Beyond the dip position, 
the height of the computed shoulder is always above the data and, thus, very likely to 
exceed also the LHC data. In comparison with other works, the height of our shoulder 
is similar to the one of [138] but almost one order of magnitude above the one of [57] . 

Considering Figs. 7.5 and 7.6 more quantitatively in the small-|t| region, one can 
use the well known parametrization of the differential elastic cross section in terms of 
the slope parameter B(s) and the curvature C(s) 



Using -B(s) from the preceding section and assuming for the moment C(s) = 0, one 
achieves a good description at small momentum transfers and energies, which is con- 
sistent with the approximate Gaussian shape of J pp (s, at small energies shown in 



da ei /dt(s,t) = da el /dt( 



s,t = 0) exp [B(s)t + C(s)t 2 } 



(7.14) 
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Fig. 6.1. The dip, of course, is generated by the deviation from the Gaussian shape at 
small impact parameters. According to (7.14), the shrinking of the diffraction peak with 
increasing energy simply reflects the increasing interaction radius described by B(s). 

For small energies \/s, our model reproduces the experimentally observed change 
in the slope at \t\ ~ 0.25 GeV 2 [157] that is characterized by a positive curvature. For 
LHC, we find clearly a negative value for the curvature in agreement with [138] but in 
contrast to [57]. The change of sign in the curvature reflects the transition of J(s, \b±\) 
from the approximate Gaussian shape at low energies to the approximate step-function 
shape (7.4) at high energies. 

Important for the good agreement with the data is the longitudinal quark momen- 
tum distribution in the proton. Besides the slope parameter, which characterizes the 
diffraction peak, also the dip position is very sensitive to the distribution width Az p , 
i.e. with Az p = the dip position appears at more than 10% lower values of \t\. In the 
earlier SVM approach [57], the reproduction of the correct dip position was possible 
without the z-dependence of the hadronic wave functions but a deviation from the data 
in the low-|t| region had to be accepted. In this low-|t| region, we achieved a definite 
improvement with the new correlation functions (2.50) and the minimal surfaces used 
in our model. 

The differential elastic cross section computed for ir^p and K ± p reactions has the 
same behavior as the one for pp (pp) reactions. The only difference comes from the 
different z-distribution widths, Az^ and Azk, and the smaller extension parameters, S n 
and Sk, which shift the dip position to higher values of \t\. This is illustrated in Fig. 7.7, 
where the model results (solid line) for the 7r ± p and K ± p differential elastic cross section 
are shown at yfs = 19.5 GeV as a function of \t\ in comparison to experimental data 
(closed circles) from [158]. The deviations from the data towards large \t\ leave room 
for odderon and reggeon contributions. Indeed, with a finite diquark size in the proton, 
an odderon occurs that improves the description of the data at large values of \t\ [159]. 

7.5 The Elastic Cross Section a el , a el /(T tot and <r tot /B 

The elastic cross section obtained by integrating the differential elastic cross section 

a " (s) = L dt nr (s - t)= L df iw |T(s ' t)|2 (7 ' 15) 

reduces for our purely imaginary T-matrix (5.15) to 

a el (s) = Jd 2 b L \J(s,\b L \)\ 2 . (7.16) 

Due to the squaring, it exhibits the saturation of J(s, \ b±\) at the black disc limit more 
clearly than o~ tot (s). Even more transparent is the saturation in the following ratios 
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Figure 7.7: The differential elastic cross section da el /dt(s,t) as a function of |t| for n^p and 
K^p reactions at a cm. energy of y/s — 19.5 GeV. The model results (solid line) are compared 
with the data (closed circles with error bars) from [158]. 



given here for a purely imaginary T-matrix 



Jd 2 b ± \J(s,\b±\)\ 2 
v tot{ 2fd% ± J(s,\b ± \) 



a tot (2jd 2 b L J(s,\b ± \) 

%r(s) = => ^- , (7-18) 

B Jd 2 b ± \b ± \ 2 J(s,\b L \) 

which are directly sensitive to the opacity of the particles. This sensitivity can be 
illustrated within the approximation 

T{s, t) = is <r tot {s) exp[5(s)t/2] (7.19) 

that leads to the differential cross section (7.14) with C(s) = 0, i.e. an exponential 
decrease over \t\ with a slope B(s). As the purely imaginary T-matrix element (7.19) is 
equivalent to 

J(s, \b±\) = {a tot /ATxB) exp[-|6 ± | 2 /2 J B] = (4<r el /0 exp[-|6 ± | 2 /2 J B] , (7.20) 
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one finds that the above ratios are a direct measure for the opacity at zero impact 
parameter 

J( a , \b ± \ = 0) = (a tot /4TrB) = (4a el /a tot ) . (7.21) 

For a general purely imaginary T-matrix, T(s,t) = i s a tot g{\t\) with an arbitrary real- 
valued function g(|i|), J(s, \b±\ = 0) is given by a el /a tot times a pure number which 
depends on the shape of g(\t\). 

We compute the elastic cross section a el and the ratios a el /a tot and a tot /B in our 
model without taking into account reggeons. In Fig. 7.8 the results for pp and pp 
reactions (solid lines) are compared with the experimental data (open and closed circles). 
The data for the elastic cross section are taken from [10] and the data for a tot and B 
from the references given in previous sections. For pp {pp) scattering, we indicate 
explicitly the prediction for LHC at yfs = 14TeV and the onset of the black disc 
limit at yfs = 10 6 GeV. The model results for up and Kp reactions are presented by the 
dashed and dotted line, respectively. For the ratios, the asymptotic limits are indicated: 
Since the maximum opacity or black disc limit governs the yfs — > oo behavior, a el /a tot 
(a tot /B) cannot exceed 0.5 (8ir) in hadron-hadron scattering. 

In the investigation of pp (pp) scattering, our theoretical curves successfully confront 
the experimental data for the elastic cross section and both ratios. At low energies, the 
data are underestimated as reggeon contributions are not taken into account. Again, the 
agreement is comparable to the one achieved in [138] and better than in the approach 
of [57], where a el comes out too small due to an underestimation of da el /dt in the low-|t| 
region. 

Concerning the energy dependence, a el shows a similar behavior as a tot but with 
a more pronounced flattening around yfs > 10 6 GeV. This flattening is even stronger 
for the ratios, drawn on a linear scale, and reflects very clearly the onset of the black 
disc limit. As expected from the simple approximation (7.21), a el /a tot and a tot /B show 
a similar functional dependence on yfs. At the highest energy shown, yfs = 10 8 GeV, 
both ratios are still below the indicated asymptotic limits, which reflects that J(s, \ b^_\) 
still deviates from the step-function shape (7.4). The ratios a el /a tot and a tot /B reach 
their upper limits 0.5 and 8ir, respectively, at asymptotic energies, yfs — > oo, where the 
hadrons become infinitely large, completely black discs. 

Comparing the pp (pp) results with the ones for 7rp and Kp, one finds that the 
results for a tot / B converge at high energies as shown in Fig. 7.8. This follows from the 
identical normalizations of the hadron wave functions that lead to an identical black 
disc limit for hadron-hadron reactions. 
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Figure 7.8: The elastic cross section a el and the ratios a el /a tot and a tot /B as a function of the 
cm. energy y/s. The model results for pp (pp), irp, and Kp scattering are represented by the 
solid, dashed and dotted lines, respectively. The experimental data for the pp and pp reactions 
are indicated by the open and closed circles, respectively. The data for the elastic cross section 
are taken from [10] and the data for a tot and B from the references given in previous sections. 
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Chapter 8 



Conclusion 



We have developed a loop-loop correlation model (LLCM) [22, 23] in which the QCD 
vacuum is described by perturbative gluon exchange and the non-perturbative stochastic 
vacuum model (SVM) [27]. This combination leads to a static quark-antiquark potential 
with color Coulomb behavior for small and confining linear rise for large source separa- 
tions in good agreement with lattice QCD results. We have computed in the LLCM the 
vacuum expectation value of one Wegner- Wilson loop, (Wj.[C]), and the correlation of 
two Wegner- Wilson loops, (W ri [C\] W r2 [C2]), for arbitrary loop geometries and general 
representations ru\ of SU(N C ). Specifying the loop geometries, these results allow us 
to compute the static quark-antiquark potential, the glueball mass, the chromo-field 
distributions of static color dipoles, the QCD van der Waals potential between two 
static color dipoles, and the 5-matrix element for high-energy dipole-dipole scattering. 
Together with a phenomenological universal energy dependence and hadron and photon 
wave functions, the latter provides the basis for a unified description of hadron-hadron, 
photon-hadron, and photon-photon reactions within the functional integral approach to 
high-energy scattering. 

We have applied the LLCM to compute the chromo-electric fields generated by a 
static color dipole in the fundamental and adjoint representation of SU(N C ). The for- 
mation of a confining color flux tube is described by the non-perturbative SVM correla- 
tions [28] and the color Coulomb field is obtained from perturbative gluon exchange. We 
have found Casimir scaling for both the perturbative and non-perturbative contributions 
to the chromo-electric fields in agreement with recent lattice data [84]. String breaking 
is neither described for sources in the fundamental representation nor for sources in the 
adjoint representation which indicates that in our approach not only dynamical fermions 
(quenched approximation) are missing but also some gluon dynamics. Transverse and 
longitudinal energy density profiles have been provided: For small dipoles, R = 0.1 fm, 
perturbative physics dominates and non-perturbative correlations are negligible. For 
large dipoles, R > 1 fm, the non-perturbative confining string dominates the chromo- 
electric fields between the color sources. The transition from perturbative to string 
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behavior takes place at source separations of about 0.5 fm in agreement with the recent 
results of Liischer and Weisz [40]. The root mean squared radius R ms of the confining 
string and the energy density in the center of a fundamental SU (3) dipole e^(X = 0) are 
governed completely by non-perturbative physics for large R and saturate as R increases 
at i?^7°° « 0.55 fm and e§^°°(X = 0)wl GeV/fm 3 . 

We have presented the low-energy theorems [42-44], known in lattice QCD as 
Michael sum rules [41], in their complete form in continuum theory taking into account 
the contributions found in [29,42] that are missing in the original formulation [41]. We 
have used the complete theorems to compare the energy and action stored in the con- 
fining string with the confining part of the static quark-antiquark potential. The com- 
parison shows consistency of the model results and indicates that the non-perturbative 
SVM component is working at the renormalization scale at which /3(g) jg = —2 and 
a s = 0.81. Earlier SVM investigations along these lines have found a different value of 
a s = 0.57 with the pyramid mantle choice for the surface [28, 29] but were incomplete 
since only the contribution from the traceless part of the energy-momentum tensor has 
been considered in the energy sum. 

A Euclidean approach to high-energy dipole-dipole scattering has been established 
by generalizing Meggiolaro's analytic continuation [45] from parton-parton scattering to 
gauge-invariant dipole-dipole scattering. The generalized analytic continuation allows 
us to derive S-matrix elements for high-energy reactions from configurations of Wegner- 
Wilson loops in Euclidean space-time with Euclidean functional integrals. It thus shows 
how one can access high-energy reactions directly in lattice QCD. First attempts in this 
direction have already been carried out but only very few signals could be extracted, 
while most of the data was dominated by noise [49]. We have applied this approach to 
compute in the Euclidean LLCM the scattering of dipoles at high-energy. The result 
derived in the Minkowskian version of the LLCM [23] has exactly been recovered includ- 
ing the well-known two-gluon exchange contribution to dipole-dipole scattering [99, 100]. 
This confirms the analytic continuation of the gluon field strength correlator used in all 
earlier applications of the SVM to high-energy scattering [30-32, 50-61]. 

The S-matrix element obtained in our approach allows us to investigate manifesta- 
tions of the confining QCD string in high-energy reactions of photons and hadrons [24] 
but leads to energy-independent cross sections in contradiction to the experimental ob- 
servation [23]. The missing energy dependence is disappointing but not surprising since 
our approach does not describe explicit gluon radiation needed for a non-trivial energy 
dependence. 

We have introduced a phenomenological energy dependence into the S'-matrix el- 
ement Sdd that allows a unified description of hadron-hadron, photon-hadron, and 
photon-photon reactions and respects S'-matrix unitarity constraints in impact param- 
eter space [23,25,26]. Motivated by the two-pomeron model of Donnachie and Land- 
shoff [62], we have ascribed to the non-perturbative and to the perturbative component 
a weak and a strong energy dependence, respectively. The constructed T-matrix ele- 
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ment shows Regge behavior at moderately high energies and describes multiple gluonic 
interactions important to respect unitarity constraints in impact parameter space at 
ultra-high energies. However, for a more fundamental understanding of hadronic high- 
energy reactions in our model, one faces the highly ambitious task to implement gluon 
radiation and quantum evolution explicitly. 

The model parameters have been adjusted to reproduce a wealth of experimental 
data (including cosmic ray data) for total, differential, and elastic cross sections, struc- 
ture functions, and slope parameters over a large range of cm. energies. The model 
parameters that allowed a good fit to high-energy scattering data are in good agreement 
with complementary investigations: The parameters of the non-perturbative component 

- the correlation length a, the non-Abelian strength k, and the gluon condensate G2 

- are constrained by lattice QCD investigations, by the string tension a of a static 
quark-antiquark pair, and by the SVZ gluon condensate G2 essential in QCD sum rule 
investigations. For the energy dependence, the exponents of the Donnachie-Landshoff 
two-pomeron fit, e so ft and ehard, have been used as an orientation for our energy ex- 
ponents e NP and e p . Besides these parameters describing dipole-dipole scattering, the 
reaction-dependent parameters in the light-cone wave functions are also consistent with 
other approaches: In the hadron wave functions, the transverse extension parameters Sh 
are in good agreement with the corresponding electromagnetic radii [61] and the width 
of the longitudinal quark momentum distributions Azh has been computed from [106]. 
In the photon wave function, the running quark masses, which coincide with the current 
quark masses for large Q 2 and the constituent quark masses for small Q 2 [52], have been 
chosen in agreement with sum rule derivations. 

Having adjusted the model parameters, we have studied S*-matrix unitarity limits 
of the scattering amplitudes in impact parameter space. On the basis of dipole-dipole 
scattering, we have shown explicitly that our model respects unitarity constraints in 
impact parameter space and, in particular, the black disc limit at ultra-high energies. 
The profile functions have been calculated for proton-proton and longitudinal photon- 
proton scattering. They show very clearly that the interacting particles become blacker 
and larger with increasing cm. energy. At ultra-high energies, the opacity saturates 
at the black disc limit while the transverse expansion of the scattered particles contin- 
ues. Increasing the photon virtuality Q 2 in longitudinal photon-proton scattering, the 
maximum opacity increases and also the energy at which it is reached for zero impact 
parameter. 

We have related the impact parameter dependent gluon distribution xG(x, Q 2 , \b±\) 
to the profile function for longitudinal photon-proton scattering and found low-x satu- 
ration of xG(x,Q 2 , \ b±\) as a manifestation of S'-matrix unitarity. In accordance with 
the profile function, xG(x,Q 2 , \b±\) decreases from the center towards the periphery of 
the proton. With increasing photon virtuality Q 2 , the increase of xG(x,Q 2 ,\b±\ = 0) 
becomes stronger towards small x and the saturation value of xG(x,Q 2 , \ b±\ = 0) in- 
creases but is reached at decreasing values of x. However, the integrated gluon distri- 
bution xG(x,Q 2 ) does not saturate because of the growth of the proton radius with 
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decreasing x observed in our approach. Similar results are obtained in complementary 
approaches [124,127,128]. 

More model dependent are the specific energies at which these saturation effects 
set in. The profile function saturates the black disc limit at zero impact parameter 
for \fs > 10 6 GeV in proton-proton scattering and for \fs > 10 GeV in longitudinal 
photon-proton scattering with Q 2 > 1 GeV 2 . In both reactions, the wave function nor- 
malization determines the maximum opacity. The saturation of xG(x,Q 2 , \b±\) occurs 
in our approach for Q 2 > 1 GeV 2 at values of x < 1CT 10 , i.e. far below the HERA and 
THERA range. 

For proton-proton scattering, we have found that the rise of the total and elastic 
cross section becomes weaker for yfs > 10 6 GeV due to the onset of the black disc limit 
at = in the profile function. This saturation of the profile function becomes even 
more apparent in the ratios a el /a tot and a tot /B which are a measure of the proton 
opacity. In contrast, no saturation effect has been observed in the slope parameter 
B(s) which is a measure for the transverse expansion of the proton. Considering the 
differential elastic cross section da el /dt, the model has described the diffraction pattern 
and also the shrinkage of the diffraction peak with increasing energy in good agreement 
with experimental data at small momentum transfers \t\. Around the dip region, where 
a real part is expected to be important, deviations from the data have reflected that 
our T-matrix is purely imaginary. Our predictions for proton-proton scattering at LHC 
(y/s = 14TeV) are a total cross section of £7*°* = 114.2 mb in good agreement with 
cosmic ray data and a differential elastic cross section da el jdt with a slope parameter 
of B = 21.26 GeV -2 , a negative curvature C < 0, and a dip at \t\ « 0.35 GeV 2 . 

For pion-proton and kaon-proton scattering, results analogous to proton-proton scat- 
tering have been obtained but with a slightly stronger rise observed in the total cross 
section. This has been traced back to the smaller size of pions and kaons in compar- 
ison to protons, S p = 0.86 fm > S w = 0.607 fm > Sk = 0.55 fm, and the perturbative 
component becoming increasingly important with decreasing dipole sizes involved. Fur- 
thermore, a weak convergence of the ratios cr^/al^ and cr^/a^ towards unity is 
predicted as the energy increases. The smaller size of the pion and kaon is reflected in 
the differential elastic cross sections da el /dt by the dip shifted towards larger values of 
\t\. 

For photon-proton and photon-photon reactions, an even stronger rise of the total 
cross section has been observed with increasing energy. As illustrated in the proton 
structure function F2(x,Q 2 ), this rise becomes steeper with increasing photon virtuality 
Q 2 . Again, we have traced back the strong energy boost to the growing importance of the 
perturbative component with decreasing dipole size. Besides some deviations from the 
experimental data with increasing Q 2 , our model has described a l °* p (s, Q 2 ) successfully 
in the low-Q 2 region where the running quark masses become constituent quark masses. 

Finally, let us emphasize that the presented Euclidean approach to high-energy scat- 
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tering is independent of our specific loop-loop correlation model. It indeed makes any 
method limited to a Euclidean formulation of the theory applicable for investigations 
of high-energy reactions. Here encouraging new results have been obtained with in- 
stantons [48] and within the AdS/CFT correspondence [46] and it will be interesting 
to see precise results from the lattice. A promising complementary Euclidean approach 
has been proposed in [160] where the structure functions of deep inelastic scattering 
at small Bjorken x are related to an effective Euclidean field theory. Here one hopes 
that the limit x — > corresponds to critical behavior in the effective theory. The aim is 
again to provide a framework in which structure functions can be calculated from first 
principles using genuine non-perturbative methods such as lattice computations. In an- 
other recent attempt, the energy dependence of the proton structure function has been 
related successfully to critical properties of an effective near light-cone Hamiltonian in 
a non-perturbative lattice approach [161]. It will be interesting to see further develop- 
ments along these lines aiming at an understanding of hadronic high-energy scattering 
from the QCD Lagrangian. 
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Appendix A 



Loop and Minimal Surface 
Parametrizations 



A rectangular loop Ci with "spatial" extension Ri and "temporal" extension Tj placed in 
four-dimensional Euclidean space — as shown in Fig. A.l — has the following parameter 
representation 



with 



where 
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The "center" of the loop Ci is given by X$i. The parameters Zi, Ri, 9i, fa, and Qi are 
defined in Fig. A.l that illustrates (a) the spatial arrangement of a color dipole and 
(b) its world-line in Euclidean "longitudinal" space. The tilting angle Qi / is the 
central quantity in the analytic continuation presented in Chap. 4. Moreover, Q\ = tt/2 
together with ©2 = allows us to compute conveniently the chromo-magnetic field 
distributions in Appendix B. 
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The minimal surface Si is the planar surface bounded by the loop Cj = dSi given 
in (A.l). It can be parametrized as follows 

Si = $yXi(ui,Vi) = Xoi - (1 - Zi)ri + Viti + UiU, m £ [-Ti,Ti], v { G [0, 1]| (A.7) 

with rj and ti given in (A. 6). The corresponding infinitesimal surface element reads 





Figure A.l: (a) Spatial arrangement of a color dipole and (b) its world- line in Euclidean 
"longitudinal" space given by the rectangular loop Ci that defines the minimal surface Si with 
dSi = Ci . The minimal surface is represented by the shaded area. In our model, it is interpreted 
as the world-sheet of the QCD string that confines the quark and antiquark in the dipole. 
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Appendix B 



%-Computations with Minimal 
Surfaces 



The quantities considered in the main text are computed from the VEV of one loop 
(W[C]) and the loop-loop correlation function (W[Ci]W[C2]). Using the Gaussian ap- 
proximation in the gluon field strengths, both are expressed in terms of xSjSj -functions 
(2.15) and (2.27) as shown in Sees. 2.1 and 2.2. These x-functions are central quantities 
since here the ansatz of the gauge-invariant bilocal gluon field strength correlator and the 
surface choice enter the model. In this Appendix these functions are computed explic- 
itly for minimal surfaces (A. 7) and the F^pu - ansatz given in (2.42), (2.44), and (2.48). 
Note that the contributions from the infinitesimally thin tube — which allows us to 
compare the gluon field strengths in surface S\ with the gluon field strength in surface 
62 — cancel mutually. 

Depending on the geometries and the relative arrangement of the loops, the \- 
functions determine the physical quantities investigated within the LLCM such as the 
static qq potential (3.2), the chromo-field distributions of a color dipole (3.21), and the 
S'-matrix element for elastic dipole-dipole scattering (4.12). 

We compute separately the three components g , X5iS™ C > an d Xs^sl ^ or g enera l 
loop arrangements from which the considered quantities are obtained as special cases. 
Without loss of generality, the center of the loop C2 is placed at the origin of the 
coordinate system, A02 = (0,0,0,0). Moreover, C2 is kept untilted, 62 = 0, and 
:= Gi is used to simplify notation. We limit our general computation to loops with 
r i,2 = (rL,2±) 0, 0) = #1,2 = vr/2 and transverse "impact parameters" 6 = Xqi — A02 = 
Aqi = (61,62,0,0) = (6i,0,0) which allows to compute all of the considered quantities. 
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X^^-Computation 



Starting with the definition 

J2 



XSist ■■= T f da^X,) [ da pa (X 2 )F^(Z = X 1 -X 2 ) (B.l) 

= 12 (jv^ 2 - 1) j s da ^ Xl ^ j s da p*( x i) (<VA<r - <WM D(Z 2 ) , 



one exploits the anti-symmetry of the surface elements, da^ v = —da u ^, and applies the 
surface parametrization (A. 7) with the corresponding surface elements (A. 8) to obtain 



7T 2 G 2 K 

ml - r 



(/ ] r 2 ) I d Vl I dv 2 I 'dm [ ' du 2 D{Z 2 ) (B.2) 

J JO J-Ti J-T 2 



with 



Z — Xi — Xo 



( b± (1 zi) r 1± + v 1 f 1± + (1 - z 2 ) f 2± - v 2 f 2± \ 
— u\ sin 
U\ cos — u 2 



(B.3) 



where the identities t± ■ r 2 = r± ■ t 2 = and t± ■ t 2 = cos 0, evident from (A. 6) with the 
mentioned specification of the loop geometries, have been used. In the limit T 2 — > oo, 
the u 2 integration can be performed 



lim [ 1 du 2 D(Z 2 ) 



d 4 k r T2 

^-,D{K 2 ) lim / du 2 e* KZ 



-I 



d 4 K ~ 



(2vr) 4 
d 3 K 
(2vr) 3 



D(K 2 )2ttS(K 4 ) exp 



iK±Z± + iK^u\ sin + iK^u\ cos 



D^ 3 \K 2 )exp 



iK±Z± + 1K3U1 sin 



= D^(Z 2 ) 



(B.4) 



which leads to 



lim X ^sl = cos© jff 2 " (fi± ■ r 2± ) f'dv, fdv 2 f^d^D^^ 2 ) . (B.5) 



T 2 

Taking in addition the limit T\ — > 00, the u\ integration can be performed as well 

rTi 



lim / du\e l 



K 3 u lsi nB = I 27^3 sin©) 
1 limr^oo 2Ti 



271^(^3 ) 

sine| 



for sin / 
for sin = 



(B.6) 



With Ti = T 2 = T/2 -► 00, one obtains for sin© / 

cos it 2 G 2 k 



T ^x SlS2 - |gin@| 3{N2 



'~~~T\ (n± ■ r 2 ±) j 1 dv 1 f 1 dv 2 D^(Zl) (B.7) 
~ l ) Jo Jo 
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and for sin 6 = 

lim j$& = lim Tcose j'f 2K (r 1± -r 2± ) fdv, fdv 2 D^(Z 2 ) . (B.8) 

Evidently, (B.7) is the result given in (4.14) and (4.17) that describes the confining 
contribution to the dipole-dipole scattering matrix element Sdd- 

From (B.8), one obtains the confining contribution to the static color dipole potential 
for S\ = S 2 = S which implies T\ = T 2 = T/2, = 0, Z\ = z 2 , r\ = r 2 = r, and 
H_L • r 2 ± = r 2 = R 2 so that 

- ^ T wf^ R2 I/ p(l - p)D ' 3,( ^ u (B ' 9) 

which leads directly to (3.7). 

From (B.5) the confining contribution to the chromo-field distributions AG 2 tj3 (X) 
can be computed conveniently. Equation (B.5) reads for Si = Sp, T\ = Rp/2 and 
Ri = R P , and S 2 = S w , T 2 = T/2 and R 2 = R 

lim X N s P J w = cos 6 J^ 2K (f 1± • r 2± ) / dv, / <fo 2 / d Ul D^(Z 2 ) (B.10) 
with 

2 = X 1 -X 2 =( b-L-( l - z i)ri_L+v 1 r 1± + (l-z 2 )r 2± -v 2 r 2± \ 
\—ui sin y 

The confining non-perturbative contribution to the chromo-magnetic fields vanishes as 
it is obtained for plaquettes with O = ir/2. The corresponding contribution to the 
chromo-electric fields can be computed with = as follows: Due to R\ = R p — > 0, 
the u\ and v\ integrations in (B.10) can be performed with the mean value theorem. 
Keeping only terms up to 0(R 2 ), the confining non-perturbative contribution to the 
chromo-field distributions AG 2 ij3 (X) is obtained as given in (3.30). 



X^^ c -Computation 



We start again with the definition 

J2 



NPnc 
X-S1S2 



J s da^X,) da pa {X 2 ) F$%?(Z = X, - X 2 ) 
K ' C hi ' / da^(Xi) [ da pa (X 2 ) (B.12) 

l ) JS! JS 2 



12{N 2 



1 

x- 
2 



d 



d 

(Z a 5 pp - ZpSpa) + (Z p 5 V( j - Z a b vp ) D 1 (Z Z ) 
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and use the anti-symmetry of both surface elements to obtain 



NPnc 
AS1S2 



TT 2 G 2 (l- 


"«) 


6(N* - 


1) 


7T 2 G 2 (1- 


"«) 


3(iV 2 - 


1) ■ 


7T 2 G 2 (1 





d 

da^ (Xi ) / da pa (X 2 ) Z CT <5 MP Di 



Si 



da pu {Xi] 



s 2 

/„ 

'S 2 



0£„ 
5 9 



with 



3(iV| - 1) J 8l 



d K jkz iV C r^2> 



^ js 2 



da pa (X 2 ) 



d 



dX- 



2a 



d 4 K lKZ d b ^ R2) 



(27TY 



dK 2 



(B.13) 
(B.14) 
(B.15) 

(B.16) 



As evident from (B.15), Stokes' theorem can be used to transform each of the surface 
integrals in Xsfs™ C a nne integral 



NPnc 
A.S1S2 



7T 2 G 2 (1- 


-K) 


3(iV 2 - 


1) 


7T 2 G 2 (1- 


-k) 


6(N 2 - 


1) 


7T 2 G 2 (1- 


-K) 


3(iV 2 - 


1) 



da fJ , u (X 1 ] 



d 



Si 



dZ, 



v JC 2 



dZ p (X 2 )b pp D' x (Z z ) (B.17) 



! da^Xi.) I dZ p (X 2 ) 5 up Z v D±(Z 2 ) (B.18) 

(B.19) 



dZ„(Xi) d dZ p (A 2 )<5 MP ^i(^ 
i J C2 

With the line parametrizations of C\ and C 2 given in (A.l) and the specification of the 
loop geometries mentioned at the beginning of this appendix, (B.19) becomes 



NPnc 



7T 2 G 2 (1- K ) 

3(iV 2 - 1) 



(B.20) 



Ti 



x { cos B / dui / du 2 

l-Ti J-T 2 



T 2 



D[ (Z 2 AA )- D[ (Zi c )- D[ {Z 2 CA )+ D[ (Z 2 CC 



V (rrl 



I trzl 



+ (r 1± -r 2± )/l 1 /^ 2 bU^ B )-^U^D)-^i(^)+^U^DD) 

JO JO L 

where the following shorthand notation is used 

Z XY := X* — X\ with X% G C% and X\ G C£" . 



(B.21) 



In the limit i?i j2 <C Ti )2 — > 00, the term oc (rij_ • r 2 j_) on the rhs of (B.20) can be 
neglected and, thus, (B.20) reduces to 



Hrn ^NPnc_ g TT 2 G 2 (l -K 



lim / dm lim / du 2 (B.22) 

1) Ti^oo7_ Tl T 2 ^ooJ_ T2 



D[(Z 2 AA ) - D[(Zi c ) - D[(Z£ A ) + £>i(Z^ c 



/ / ry 2 
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Here the integrations over u\ and u 2 can be performed analytically proceeding analo- 
gously to (B.4) and (B.6). With T x = T 2 = T/2 -»• oo, one obtains for sin 9 / 



lim v NPnc 

lim x 5l& 

1 — »oo 



cosG ir 2 G 2 (l-n) 



|sin9| 3(iV 2 -l) 
-^1 \ /j AA±)~ u l \ Zj AC±)~ u l \ Zj CA±) + u l \ A CC±) 



(B.23) 



and for sin 9 = 



lim ^ NPnc 

lim XSl s 2 

1 — >oo 



7T 2 G 2 (1 - K) 



lim T cos 6 ■ . , . , 

T^oo 3(N 2 - 1) 



(B.24) 



With the following identities 

Zaal = n q -r 2q , Z AC± = r lq -f 2 q , Z CA± = r lq —r 2q , Z C c_L = r lq —r 2q - 



(B.25) 



one sees immediately that (B.23) is the result given in (4.14) and (4.16) that describes 
the non-confining non-perturbative contribution to the dipole-dipole scattering matrix 
element Sdd- 

From (B.24) one obtains the non-confining contribution to the static potential for 
S 1 = S 2 = S, i.e., T 1= T 2 = T/2, 9 = 0, n = r 2 = r, 



lim Xss 

1 -^oo 



NPnc 



- lim T 



tt 2 G 2 (1-k) 



X 



r->oc 3(N 2 - 1) 
DT\Z\ A ) - Df\z\ c ) - Df\zl A ) + Df\zl c ) 



(B.26) 



which contributes to the self-energy of the color sources with 

r 2 



lim -v-^Pnc 
„ ilm XSS self 
1 — >oo 



Km T - i> 

T^oo 3(iV 2 _ 1) 

2^ 2 G 2 (1-k) 

hm 1 



T^oo 3{N 2 - 1 

and to the potential energy between the color sources with 

^ 2 G 2 (1 - « 



dT\z\ a ) + Df\z 2 c) 

dT\z\ a ) 



(B.27) 



lim yf;5 = lim T 

T^oo ASSpot T^oo 6(iV|-l) 



= lim T 



tt 2 G 2 (1-k) 



dT\z\ c) + d^\zI a ) 

D?\z\ c ) . 



(B.28) 



r^oo 3(iV 2 - 1) 

The latter gives the non-confining contribution to the static potential (3.6). 

From (B.18) the non-confining non-perturbative contribution to the chromo-electric 

T af3\ 



fields (AG 2 t[3 (X) with a(5 = z4 = 4z) can be computed most conveniently with zero 
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plaquette tilting angle = 0. The corresponding contribution to the chromo-magnetic 
fields (AG^g(X) with a(3 = ij = ji) is obtained for plaquette tilting angle = it/2 
and thus vanishes which can be seen most directly from the surface integrals (B.13). 
Now, we set = to compute the contribution to the chromo-electric fields: Using 
the surface Si = Sp and loop C2 = dSyy parametrizations, (A. 7) and (A.l), with our 
specification of the loop geometries, one obtains from (B.18) 



NPnc 
X-SpSw 




R P /2 rl 
Rp/2 JO 



(fii. ■ Z 1A± ) D X {Z{ A ) - (f 1± ■ Z 1C± ) D^Zfc) 



(B.29) 



dV2 



(fi± ■ Z 1B± ) D±{Zl B ) - (r 1± ■ Z 1D± ) D x {Z\ Di 



-T/2 

- (n± ■ f 2 ±) / a 
Jo 

with Ti = Rp/2, Ri = R P , T 2 = T/2, R 2 = R, and the shorthand notation 

Z 1X := Xi - X% with X\ 6 Si = S P and X? € Cg = dS^ . (B.30) 
In the limit R <C T — > 00, the term oc (fi± ■ r*2j_) on the rhs of (B.29) can be neglected 

-Rp/2 



lim XSpS 

1 — >oo 



NPnc 
pSw 



— -f / dm / dvi lim / du 2 (B.31) 

- 1) J-Rp/2 Jo T -"x>J-T/2 



3( N c ~ L ) J-Rp/2 JO *->°oJ-T/2 

{ri± ■ Zia±) D x (Zl A ) - (r 1± ■ Z 1C ±) B X {Z\ C ) 



With an infinitesimal plaquette used to measure the chromo-electric field, R\ = R p — > 0, 
the mean value theorem can be used to perform the u\ and v\ integrations in (B.31). 
Keeping only terms up to O(Rz), this leads directly to the non-confining non-per- 
turbative contribution to the chromo-field distributions AG^g(AT) as given in (3.28) 
and (3.29). 



X P -Computation 



Comparing the definition of the perturbative component 

Xis, ■■= ^j s da^(X 1 )j s da pa (X2)F^ pa (Z = X 1 -X2) 

= £ / tfcvpTi) / da pa (X 2 ) (B.32) 
4 JSi Js 2 

1 r d d 1 

x 2 ~qz~ { z <t$w ~ z p$ua) + -Qg- (Z p 5 ua - Z a 5 up ) D P (Z ) 

with the one of the non-confining non-perturbative component Xsfs™ gi ven m (B.12), 
one finds an identical structure. Thus, accounting for the different prefactors and the 
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different correlation function, the results for Xs 1 s 2 can ^ e reac ^ on? directly from the 
results for Xsisl° gi ven above: 

With T\ = T2 = T/2 — > 00 and our specification of the loop geometries, one obtains 
the result for sin / from (B.23) 

, . p cos 

lim XSl S 2 = ~ 



1- - x I sin 9 

i'(2)/^2 \_ n '( 2 )f72 \_n'(2)^2 \^-n'W(7- 



(B.33) 



X 



Dp ( Z AA±)-D P (Z AC ±)-D p {Z CAA _) + D p (Z cc± ) 



and the result for sin 6 = from (B.24) 



lim X Sl S 2 

1 — >oo 

y(3)/ 



lim T cos q 



(B.34) 



n > (3)/ 7 2 



+ < ) (4c) 



where Z^y is defined in (B.21) and Zxy± is given explicitly in (B.25). Evidently, (B.33) 
is the final result given in (4.14) and (4.15) that describes the perturbative contribution 
the dipole-dipole scattering matrix element Sdd- 

The perturbative contribution to the static potential is obtained from the expression 
corresponding to (B.26), 



lim Xss = ~ I™ T f 2 

1 -^OO 1 ~^00 



(B.35) 



Df\z AA ) - D' p ^{Z\ c ) - D'^{Z 2 CA ) + n'?\zl c 
which contributes to the self-energy of the color sources with 

Df\z\ A )+Df\zl c ) 



y(3) 



'(3), 



y(3), 



lim xfs S eif = - lim T g 2 

1 -^00 1 — >oo 



= - lim T2g< Df\z\ A ) 
and to the potential energy between the color sources with 



(B.36) 



lim Xss pot = - lim T g 2 

1 — >oo 1 -^00 



Df\z\ c )+Df\zl A ) 



= - \^T2g 2 Df\z\ c ) . 
The latter gives the perturbative contribution to the static potential (3.4). 



(B.37) 



The perturbative contribution to the chromo- magnetic fields (AG 2 t/3 (X) with a(3 = 
ij = ji) vanishes while the one to the chromo-electric fields (AG^(X) with aj3 = 
i4 = 4i) for which a plaquette with = is needed, is obtained from the expression 
corresponding to (B.31), 

r-T/2 



rRp/2 rl rT/2 

lim Xs P s w = -9 2 / du i / ^1 lim / du 2 

T ^°° J-Rp/2 JO T-"x>J-T/2 



(B.38) 



(rii. • Z 1A± ) D P {Z 2 lA ) - (f 1± ■ Z 1C± ) Dp(Z 2 c ) 
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with Z\x as defined in (B.30). To perform the u\ and v\ integrations in (B.38), again the 
mean value theorem can be used since the plaquette has infinitesimally small extensions, 
R\ = R p — > 0. Keeping only terms up to O(Rp), this leads directly to the perturbative 
contribution to the chromo-field distribution AG^JX) as given in (3.25) and (3.26). 
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